Spin-1 effective Hamiltonian with three degenerate orbitals: 
An apphcation to the case of V2O3. 
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Motivated by recent neutron and x-ray observations in V2O3, we derive the effective Hamiltonian 
in the strong coupling limit of an Hubbard model with three degenerate t2g states containing two 
electrons coupled to spin S = 1, and use it to re-examine the low-temperature ground-state prop- 
erties of this compound. An axial trigonal distortion of the cubic states is also taken into account. 
Since there are no assumptions about the symmetry properties of the hopping integrals involved, 
the resulting spin-orbital Hamiltonian can be generally applied to any crystallographic configuration 
of the transition metal ion giving rise to degenerate t2g orbitals. 

Specializing to the case of V2O3 we consider the low temperature antiferromagnetic insulating 
phase. We find two variational regimes, depending on the relative size of the correlation energy of 
the vertical pairs and the in-plane interaction energy. The former favors the formation of stable 
molecules throughout the crystal, while the latter tends to break this correlated state. Using the 
appropriate variational wave functions we determine in both cases the minimizing orbital solutions 
for various spin configurations, compare their energies and draw the corresponding phase diagrams 
in the space of the relevant parameters of the problem. We find that none of the symmetry-breaking 
stable phases with the real spin structure presents an orbital ordering compatible with the magnetic 
space group indicated by very recent observations of non- reciprocal x-ray gyrotropy in V2O3. We 
do however find a compatible solution with very small excitation energy in two distinct regions of 
the phase space, which might turn into the true ground state of V2O3 due to the favorable coupling 
with the lattice. We illustrate merits and drawbacks of the various solutions and discuss them in 
relation to the present experimental evidence. 

PACS numbers: 71.27.-f-a, 71.30. +h, 75.25.+Z 



I. INTRODUCTION 



The crystal and electronic structure of V2O3 (vana- 
dium sesquioxide) has been the subject of intensive the- 
oretical rftwi experimental studies over the past three 
decades.ErEl As is well known, this compound is con- 
sidered to be the prototype of the Mott-Hubbard sys- 
tems, showing metal-insulator transitions from a para- 
magnetic metallic (PM) phase to an antiferromagnetic 
insulating (AFI) phase at low temperatures (~ l^QK) 
and from a PM phase to a paramagnetic insulating (PI) 
phase at higher temperatures (« 500 A'), due to the in- 
terplay between band formation and electron Coulomb 
correlation.B Actually it is the only known example 
among trsmsition-metal oxides to show a PM to PI 
transition.El In the paramagnetic phases the crystal has 
the corundum structure, in which the V ions are arranged 
in V-V pairs along the c-hexagonal axis and form a hon- 
eycomb lattice in the basal ah plane. Each V^+— ion has 
?)<P configuration and is surrounded by an oxygen octa- 
hedron with a small trigonal distortion which lifts the 
three-fold degeneracy of t2g orbitals into a nondegener- 
ate aig and doubly degenerate eg orbitals separated by 



the distortion energy At, with the Cg states lying lower. 
Cooling down to the AFI phase, a strongly destructive 
first order transition takes place and the system becomes 
monoclinic. At the same time a peculiar antiferromag- 
netic (AF) spin structure that breaks the original trig- 
onal symmetry of the corundum in the basal plane sets 
in. The magnetic order consists of ferromagnetic planes 
perpendicular to the monoclinic bm axis that are stacked 
antiferromagnetically. This is rather surprising, since the 
corundum lattice is known to be non frustrated. 

Different theoretical models have been proposed to ex- 

3i&¥(CNR) at- 
tempted in a series of papers a realistic description of 
the complex magnetic properties and phase diagram of 
V2O3. They realized that the pecuhar structure observed 
in the AFI phase could not be explained in terms of a 
single band Hubbard model and that the introduction of 
additional degrees of freedom in the model, in terms of 
orbital degeneracy of the atomic 3d states involved, was 
necessary to explain the experimental findings. They ar- 
gued that there was only one magnetic electron in the 
doubly degenerate eg band, the other electron being in- 
volved in a strong covalent diamagnetic bond formed by 



plain the anomalous properties of V2O3P r^n the late 
seventies Castellani, Natoli and Ranninger 



the flig states along the vertical pair. In this way a 
phase diagram showing the regions of stability of the 
various spin-orbital orderings was obtained as a func- 
tion of the in-plane hopping integrals and the J /U2 ra- 
tio, J being the intra-atomic exchange integral and U2 
the Coulomb integral between different orbitals. Implicit 
in the model was the assumption that J was not too 
strong, so that a spin S ~ 1/2 for the l/'^+— ion was 
favored instead of S* = 1. The importance of the or- 
bital degrees of freedom in the physics of V2O3 was later 
confirmed by several experiinental facts. A series of neu- 
tron scattering experimentsQiiy showed evidence that the 
wave vector of the short range antiferromagnetic corre- 
lations in both paramagnetic phases was rotated by 30 
degrees with respect to that observed in the AFI phase. 
In all three phases the atoms along the vertical pairs 
are ferromagnetically coupled; however one passes from 
a complete antiferromagnetic coupling along the three 
in-plane bonds in the high temperature paramagnetic 
phases to one ferromagnetic and two antiferromagnetic 
in the low temperature phase, with consequent breaking 
of the corundum trigonal symmetry. Moreover the mag- 
netic fluctuations in the PM and PI phases remain short 
range down to the transitiop-|temperature. Finally NMR 
studies by Takigawa et al^ in the PM phase showed 
clear evidence of the role of the orbital degrees of free- 
dom in the relaxation mechanism of the nuclear spins. 
These findings strongly point to an interpretative scheme 
in which the orbital degrees of freedom are frozen in the 
AFI phase to give rise to the peculiar AF spin structure 
and are responsible for the short range magnetic fluctu- 
ations in the high temperature phases. 

However a direct experimental evidence of the orbital 
ordering (00) and of the spin state of the V-ions was 
still lacking and was made possible only by a series of x- 
ray spectroscopies carried out with the very bright beam 
of a third generation synchrotron radiation source. For 
the first time one could subject to experimental test the 
assumptions and the predictioiis, of the CNRQ model. It 
all started when Fabrizio et alcn suggested thatJ-he ex- 
istence of an 00 of the kind suggested by CNRD could 
be revealed by means of x-ray resonant diffraction at the 
VciBadium K-edge. Following this suggestion Paolasini et 
aZ.Lj carried out such measurements. On the basis of a 
classification of forbidden reflections into pure magnetic 
and orbital, they interpreted the (111) monoclinic reflec- 
tion as evidence of an orbital ordering in V2O3, in keeping 
with the prediction of the CNR model. However in the 
same paper a measurement of the ratio between orbital 
and spin moments by non resonant magnetic scattering 
provided a value {L)/{S) = —0.3 which, together with 
the value of the magnetic moment (L) + 2{S) — 1.2 /is 
seen by neutrons,Q gives 2(5') = 1.7/1^, compatible more 
with a spin S = 1 than with a spin S = 1/2 state of the 
V atoms. This finding was pointing to an inconsistency 
of the interpretative framework. Additional evidence for 
a spin 5=1 state of the Vanadium atoms came from 
the interpretation of linear dichroism in the Vanadium 



L-edge absorption spectra|lil where evidence of a reduced 
occupation of aig orbitals was also found (25% in the PM 
phase, 20% in the PI phase, 17% in the AFI phase, in- 
stead of the 50% postulated by CNR). These are strong 
indication that in the AFI phase of V2O3 intra-atomic 
correlations prevail over band delocalization, contrary to 
the assumption made by CNR.Q 

Actually the transition between the two regimes had 
been examined in the second paper of the series CNR,eI 
where a realistic calculation was performed in the 
Hartree-Fock (HF) approximation, using a bare tight 
binding band structure for the aig and Cg orbitals in the 
Hubbard Hamiltonian. This method is very akin to the 
more modern LDA + U mcthod.li3 There it was shown 
that the value J/Ui « 0.2, where Ui is the intra-atomic 
Coulomb parameter within the same orbital, marks the 
transition between a V spin 5=1/2 regime described 
above and a spin 5 = 1 regime, where a stable self- 
consistent solution was found with the real (i.e., ob- 
served) spin structure (RS), no orbital ordering, 1.5 elec- 
trons in the eg and 0.5 in the aig bands aligned by intra- 
atomic exchange, with a spin moment of 1.7/is. This 
solution was metallic but was stabilized by the mono- 
clinic distortion with the concomitant opening of a gap 
for convenient, reasonable values of U and J. At the time 
this possibility was discarded in favor of the spin 5=1/2 
solution showing an orbital ordering, mainly ouj-the basis 
of general considerations of the phase diagram.a 

More pfficently Ezhov et alS3 using the LDA +U 
method,t3 proposed substantially the same solution with 
the same stabilization mechanism (with no surprise since 
the two methods are conceptually identical and the range 
of parameters and the initial bare density of states rather 
similar). Relying on an estimate of the effective electron- 
electron interaction parameters performed by Solovyev et 
aLili on the basis of the LDA + U method in perovskite 
LaVOa, they noted that J w 1 eV, not screened in the 
solid, and C/2 ~ 3 eV. Moreover their LDA calculations 
indicate that the aig — Cg splitting due to the trigonal 
field is At w 0.4 eV with the Cg band center below that 
of the aig band. They then argued that with these values 
of J and At the IsgCg) configuration for the AFI is more 
favorable than the |aige^) configuration, implying that 
the AFI ground state is not degenerate (no 00) with 
spin 5 = 1 on each V atom. 

However the fact that the spin 5=1 solution seems 
to point to a lack of orbital ordering must be an artifact 
of the Hartrec-Fock approximation, implicit in the LDA 
+ U scheme, together with the high J value, since an ex- 
amination of the electronic states of the F^+— ion leads 
again to an orbital degeneracy. Indeed out of the three 
one electron states |eg), jCg) and |aig) in octahedral sym- 
metry one can form three degenerate two electron states: 
Kel) = |0), Jaigei) = |-1), Wige^ = |1) (they con- 
stitute the spin and orbital triplet ground state of a two 
electron system in a strong cubic crystal fieldE3). In the 
presence of a trigonal distortion At > 0, the singlet state 



|0) would lie lowest, followed by the doublet | + 1), | — !)■ 
Putting the two atoms belonging to a vertical pair in the 
sairifi lejcg) state to gain A*, as suggested by Ezhov et 
al.f^ would loose the much bigger gain coming from the 
allowed hopping processes in a configuration in which one 
atom is in the state le^eg) and the other in \aige 
vice- versa. 
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.situation was realized by Mila et al. in their 
pape]£j where they took up an old suggestion by J. 
AUenEij that "the magnetic and optical properties of all 
the phases of V2O3 show a loss of y^+-ion identity". 
These findings, together witli|th|e results by inelastic neu- 
tron scattering quoted aboveH'li3 indicate that the vertical 
bond is quite stable and coupled to total spin Stot — 2 
with the non polar part of the wavefunction given by 
|V'±i) = i\0)a\ ± 1)& + I ± l)a\0)b)/V2 (where a and b 
indicate the two V centers). This state is clearly dou- 
bly degenerate, due to the freedom in the choice of the 
two degenerate states | ± !)■ As a consequence Mila et 
al. proposed a simple spin-orbital Hamiltonian in which 
the vertical pairs are described by a spin Stot = 2 and a 
pseudospin r = 1/2 for the orbital degeneracy, this last 
being lifted by the in-plane interaction of the pairs. The 
result was a phase diagram in which the spin structures 
with three in-plane AF bonds, with two AF and one fer- 
romagnetic (F) bond, with one AF and two F bonds and 
finally with three F bonds (respectively called G, C, A, 
and F in their paper) follow each other as being the most 
stable solutions as a functian. of the increasing parameter 
J/Ui (Fig. 2 of Mila et alJtB). Moreover spin structures 
C and A present a ferro-orbital configuration, in the sense 
that the vertical pairs are found in the same orbital con- 
figuration throughout the crystal, |V'-i) for C and IV'+i) 
for A. Structure C is indeed the one actually realized in 
the AFI phase of V2O3. This model seems to reconcile 
the existence of orbital degeneracy with the spin S — 1 
state of the V-ions. However, as will be shown below, 
there still remain problems. The stability region of the 
C phase in the parameter space is very small, like its 
stabilization energy and the percentage of occupancy of 
the \aig) with respect to the \eg) states in the molecular 
state |V'±i) of the AFI phase is 25%, t[0,be compared with 
the value of 17% found by Park et a/.til Furthermore, the 
magnetic group corresponding to the ferro-orbital order- 
ing of the C phase is not compatible with that suggested 
by the very recent observation of non-reciprocal x-ray gy- 
rotropy by Goulon et aln3 Finally, this solution cannot 
give rise to mixed orbital and magnetic reflections, as 
claimed by Mila et ali3, but only to pure magiaetic, so 
that it is also in conflict with Paolasini et alt3'c3 data, 
as will be discussed in section VII. 

All these considerations led us to re-examine the micro- 
scopic Hubbard Hamiltonian used to describe the ground 
state properties of V2O3 and to study the strong coupling 
limit of this Hamiltonian with three bands and two elec- 
trons per site couolcd to spin S — 1, along the patterns 
developed in CNRQ for spin S — 1/2, in the hope that 



starting from the fundamental Hamiltonian one could 
cure the probiems still present in the model proposed 
by Mila et altB This model is indeed based on interact- 
ing vertical V-pairs coupled to spin Stot = 2 and it is not 
clear at present what is the relation of their spin-orbital 
Hamiltonian to that derived from the atomic limit. 

The paper is organized as follows. In sections II and 
HI we introduce the three-band Hubbard Hamiltonian 
and perform a second order perturbation expansion in 
the hopping parameters in order to derive an effective 
5 = 1 Hamiltonian (-ffeff) for a system with two par- 
ticles in three degenerate orbitals. Introducing a pseu- 
dospin representation for the orbital degrees of freedom, 
we rewrite this Hamiltonian in terms of spin and pseu- 
dospin operators. The resulting spin-orbital Hamilto- 
nian, whose novelty consists in the fact that both spin 
and pseudospin degrees of freedom are 3-components vec- 
tor operators {S , l . t — 1), complements those intro- 
duced for cuprates,c3lHJ for V2O3I3 and more recently for 
manganites£3~t3 Since no assumptions are made about 
the crystal symmetry for the hopping parameters, this 
Hamiltonian can be generally applied to any crystallo- 
graphic configuration of the transition metal ion giving 
rise to degenerate t2g orbitals. We also discuss some gen- 
eral facts about Hcs, prove a theorem regarding the spin 
ground state for each bond and generalize it to the quasi- 
degenerate case, where a crystal field splitting is present. 

In view of an application to V2O3, Section IV presents 
the crystallographic and magnetic structure of this com- 
pound and discusses the magnetic space group, both in 
relation to neutron measurements and to recent observai 
tion of non-reciprocal x-ray gyrotropy by Goulon et al.c3 
From these experiments we are able to determine the true 
magnetic space group of the AFI phase, which should co- 
incide with the invariance group of the broken symmetry 
phase derived from the minimization procedure of HeS 
over the appropriate variational state. This section also 
fixes the parameters of the Hubbard Hamiltonian which 
are thought to be appropriate to the case of V2O3. 

Preliminary to the actual minimization for the entire 
crystal. Section V is devoted to an in depth study of the 
energetics of the vertical pairs, both from the point of 
view of the the original Hubbard Hamiltonian and the 
effective spin-orbital Hamiltonian. Apart from checking 
that they give the same result in the atomic limit, we 
also establish the limit of validity of Heft- Finally, this 
study will be used to determine the regions of stability of 
the various competing molecular states in the parameter 
space described by the hopping integrals, J/U2 and the 
trigonal distortion Aj, in order to assess the best varia- 
tional wave function for the entire crystal. 

We find two variational regimes and they are consid- 
ered in Section VI: in the first one, the correlation en- 
ergy of the vertical pairs is big compared to the in-plane 
interaction energy and this favors the formation of sta- 
ble vertical molecules throughout the crystal (molecu- 
lar regime). In the other situation, the bigger in-plane 



interaction energy favors uncorrelated atomic sites, due 
to a larger variational Hilbert space, thus breaking the 
molecules (atomic regime). In this Section we specialize 
the effective Hamiltonian derived in Section III to the 
case of V2O3 in the AFI phase, to determine its orbital 
and spin ground state phase diagram in the two regimes. 
In the molecular regime, besides the.spin-orbital configu- 
rations already found by Mila et alS3 (phases A, C, G, F), 
a new stable phase (RS') appears with the same real spin 
structure as the C phase and an anti-ferro in-plane orbital 
ordering, but the magnetic groups associated to the C 
and RS' phases are incompatible with the one suggested, 
by the non-reciprocal x-ray gyrotropy measurements.Ell 
We do however find a compatible solution in the same 
region of the phase space as the RS' solution with an ex- 
citation energy of less than 1 meV. In the atomic regime 
we still find a solution with the real spin structure, how- 
ever with an orbital ordei: again incompatible with the 
findings by Goulon et a/.EJ 

Finally Section VII reviews the implications of the re- 
sults obtained in the previous sections in relation to the 
present experimental evidence and provides an outlook 
on still open problems. 

Since the derivation and the form of the effective 
Hamiltonian is rather cumbersome, for convenience of 
the reader we have tried to use, whenever possible, a pic- 
torial representation of the relevant states, deferring the 
actual calculations and the final expressions to Appen- 
dices A to C. Appendix D contains useful formulas to 
calculate the average of iJoff over molecular variational 
wave functions for various spin orderings. 



II. THE MODEL 
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The meaning of these five on-site terms is the following: 

• iJg ■ describes the Coulomb repulsion between two 
electrons on the same orbital; 

(2) 

• iJg represents the repulsion between electrons 

with the same spin on different orbitals, given by 
the Coulomb minus the exchange energy; 

(3) 

• Hq, describes the Coulomb repulsion between two 
electrons with opposite spin on different orbitals; 

• Hq- represents the energy due to the jump of a pair 
of electrons with opposite spins from one orbital to 
another; 



To be as general as possible, in the present and in 
the following sections we shall ignore the specific crys- 
tallographic symmetries of V2O3 and simply deal with 
the strong coupling limit of a Hubbard Hamiltonian with 
three degenerate i2g states containing two electrons cou- 
pled to spin S = 1. This means that we shall treat 
the intra-atomic exchange J as an high energy param- 
eter with respect to the hopping integrals. The quasi- 
degenerate case, in which the three degenerate t2g levels 
are split by a trigonal field of the same order magnitude 
as the hopping term, can be easily accommodated in the 
formalism. _ 

Following CNRQ for the notations, we work with the 
same trigonal basis Cg \ eg and oig, to be referred in the 
following as orbitals I, 2 and 3. The total Hamiltonian 
can then be written as: 



H = Ho + H' 



(2.1) 



where Hq consists of the sum over the whole lattice of the 
atomic interaction terms i?o" (where j is the site-index 
which for simplicity of notation will be dropped from the 
fermion operators): 



• Hq-' is the exchange term of the process described 



in H, 



(3). 
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The kinetic (ifj) and crystal field {H'^j) terms are given 
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H't + H'^f , 













(2.3) 



where the summation is over all sites j and all possible 
spin (tj =t,i) and orbital {m,m' = 1,2,3) configura- 
tions. Moreover Ai = A2 = and A3 = A* > 0. 

Since the trigonal field splits the degeneracy of the two 
electron states by an amount comparable with the hop- 
ping integrals, in performing the atomic lirait we shall 
apply quasi degenerate perturbation theory£j whereby 
Hq is the unperturbed Hamiltonian with the reference 
energy of all three levels equal to zero. The perturbation 
term H' will then lift the spin and orbital degeneracy of 
the ground state of Hq. Our aim will be to find out a 
representation of this perturbation Hamiltonian in terms 



of the spin and orbital degrees of freedom, to describe 
the insulating phase-yf a spin-1 Mott-Hubbard system. 

As shown in CNRQ in the trigonal basis we still have: 
f/ii = U22 = U33 = Ui; f/12 = [/i3 = U23 = U2; 
J12 — J23 — J13 = J, together with the relation 



C/i = C/2 + 2 J 



(2.4) 



Consider first the zeroth-order ground state of H. The 
strong Hund's coupling J favors the triplet states with 
energy Et = U2 — J with respect to the singlet with en- 
ergy Es =1/2 + J, and to the states with two electrons on 
a single orbital, with even higher energy Ed = f/i + J = 
U2 + 3 J. Hence, in keeping with the above assumptions, 
we can construct our zeroth-order Hilbert subspace using 
only triplet states with energy Et, dropping out all the 
others. 

The total ground state of Hq can be written as a tenso- 
rial product of the atomic states over the entire crystal: 



l*n) = 



N 
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where \aj) denotes the 9-fold degenerate atomic state on 
site j with spin S = 1 and N is the total number of 
sites. Therefore |\l/o) is 9^-fold degenerate. The atomic 
subspace can be pictorially represented as follows: 
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In this notation, the circle represents a site and each 
sector in a given circle denotes an orbital, according to 
the following prescription: 




For example, |q;)i = cj^c^^ |0). 

We can now introduce the perturbation H' that par- 
tially removes the degeneracy of the \aj). Then the quasi- 
degenerate perturbation theory up to second order gives 
the following eigenvalue equation: 



^ {a\H',\(3){(3\H't\c^') 
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(2.5) 
L^f is site diagonal and H't changes the site oc- 
cupation. Here \a) and \a') are two particular states 
belonging to the 9^ degenerate ground state manifold 
and 1/3) is one of the intermediate states with one site 
singly and another site triply occupied. The first order 
term is trivial and only partially lifts the degeneracy of 
the \a) states according to their orbital population. It 
will be taken into account at the end. We shall therefore 
concentrate on the second order term. Since H^. involves 
only two sites in the hopping, the difference between the 
excited state |/3) and the zeroth-order ground state \a) 
is only in these two sites. For both sites the atomic \aj) 
is a two-electron state, while one of the atomic \f3) state 
is one-electron and the other is a three-electrons state. 
This implies that the denominator Ea — Ep, which is the 
energy difference between initial (|a)) and intermediate 
(1/3)) states for the whole crystal, is actually given only by 
the contribution of the two sites involved in the hopping 
process. Thus in the eigenvalue equation (2.5), the en- 
ergy of the ground state should be taken as the energy of 
two sites: E^ = 2Et = 2{U2 — J)- On the other hand, the 
energy Ep consists only of the three electrons-site contri- 
bution, as the one electron atom does not contribute to 
Ho- 

We consider the full multiplet structure of the interme- 
diate states, |/3a), i.e., the eigenstates of Hq with three 
electrons. They are twenty, ten of which are shown below 
(only the site with three electrons) and the other ten are 
obtained simply by reversing the spin. 
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can contribute. In this same subspace, the operator 
CjmcrXjC^^,^, is equivalent to: 
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Due to the cubic symmetry some of the intermediate 
states are degenerate, so that there are only three differ- 
ent excited levels with energies 



Efja = Ef3, — 3(C/2 - J) , 

EjSi = Ep^ = Ep^^ = Epg = Epg = 3C/2 
Ep2 = Ep^ — Epg — 3U2 + 2 J . 



(2.6) 



Accordingly, this implies that we have only three differ- 
ent values for the denominator: 



Ea ~ Ep^ — 



-(C/2-J) 

-(C/2-h2J) 

-(C/2 + 4J) 



(2.7) 



Another classification of the states can be done accord- 
ing to their total spin: it is easy to check, by summing 
three spin i, that |/3o) has 15 = |,5^ = |), l/Sy) has 
IS* = |, S'z = i) and all the others have \S = i, 5"^ = i). 
The same criteria apply to the other 10 states with op- 
posite spin, by reversing the sign of Sz- 

To proceed further, we introduce the following opera- 
tor: 
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Ea — Ep 



which can also be written more explicitly: 
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where X- ' collects the subspace spanned by the eigen- 
vectors (|/3o)j l/^y)) and those with the same orbital oc- 

(2) 

cupancy, but with opposite spin, A" ■ corresponds to the 

eigenvectors (|/5i), l/Jg), l/Ss), l/^s), l/Sg)) and xf to 
(1/^2), I/34), l/^e))] plus those with opposite spin, respec- 
tively. Explicit expressions of the operators A ■ are pre- 
sented in Appendix A. 

In the subspace spanned by the ground states \a) only 
terms of the kind 



j.nrrij.m'n' + V.^+ /.. , , 



(2.10) 



The last expression is used to reduce the number of 
fermion operators on site j. The commutators are eval- 
uated in Appendix B. 

The explicit form obtained for the effective Hamilto- 
nian H^q, expressed in terms of the fermion operators, is 
reported in Appendix A. It is composed of three terms: 

(3) 

H^fl , corresponding to those in Eq. 



H'^+H^^ 



& 

Our next task is to find a representation which allows 
to rewrite H^g in terms of two operators describing, re- 
spectively, the spin and the orbital degrees of freedom. 



III. THE SPIN-ORBITAL REPRESENTATION OF 
THE EFFECTIVE HAMILTONIAN. 



We can characterize each atomic state \aj) by two 
quantum numbers: the spin 5" and pseudospin t. The 
pseudospin operator describes the orbital occupation and 
has an algebra which is exactly analogous to that of the 
usual spin operator. Due to the 9-fold degeneracy of 
each \aj) state, we need a representation with both total 
Sj — 1 and Tj ^ 1. 

Consider the pseudospin representation, first. The or- 
bital quantization axis can be selected arbitrarily, and we 
choose the following convention: 




C3+C+IO) 



- ^+^+1 



44lo) 



|t = 1,t, = -1) 



k = i,T, = 0) 



(3.1) 



Note that this representation is valid for both spin di- 
rections, so that we can omit the spin indices, and con- 
sider only the orbital ones: c+ (m = 1, 2, 3). 

With this choice we have the following relations be- 
tween the fermion and the pseudospin operators: 



T+ = ^/2(c+c3-c+Cl), 



T = \/2(c+C2 - Cj^Cs) , 



= cjc2 - C+Ci 



(3.2) 



(3.3) 



(3.4) 



The factor V2 is necessary in order to have the correct 
commutator, [t+,t^] = 2tz, and the correct value of the 
matrix element for a spin one operator, namely: 



^ma'^ma-C^,-^C„iiicr 



i(l-52) 



'^mT'^'"T'^m"T'^™'T = "2'('^^ + 1) 

CmiC™4C+,^c„,"| = -^S^S- : |1,0) ^ |1,-1) 



|1,0)^|1,0) 

|1,0)^|1,0) 

11,1)^11,1) 



^ml'^"iiC^r^Cm"l 



,n =-hS+S,:\l,-l) 



|1,0) 



C;^|Cm|C+,j^Cm"| = \S S ^ : |1,1) ^ |1,0) 
CmTC^'tCm"! = 5'S'2S'+ : |1,0) -^ |1,1) 



(3.8) 






T±|l,r,) = V(lTT.)(2±T,)|l, r,±l) 



(3.5) 



All possible orbital transitions can be described in the 
following way with the pseudospin operators: 

C+C1C+C2 = (l + r,)(l -r,) : |1,0> ^ |1,0) 
C+C1C+C3 = -iT,(l - rj : |1, -1) ^ |1, -1) 
C+C2C+C3 = +iT,(l + r,) : |1, 1) ^ |1, 1) 
c+Cic+C3 = -^T+r, :|1,-1)^|1,0) 
C+C1C+C2 = -^r,T- : |1, 0) ^ |1, -1) 
c+C2C+C3 = -^r-r, :|1,1)^|1,0) 

C+C2C+Cl = -^T,T+:|l,0)^|l,l) 

C+C3C+C2 = +^T-r- : |1, 1) ^ |1, -1) 



(3.6) 



C3 C3C2 Cl 



T+T+ : 11,-1) 



|1,1) 



Equivalence in Eq. ( |3.6| ) between fermion and pseu- 
dospin operators must be interpreted as equality of ma- 
trix elements between corresponding states according to 
the representation (3.1). The states after colons indicate 
the only matrix element different from zero. For example, 
the fourth line means that: 



{a\ctcictc3\a') 



1 



(l,r.| -^r+r. 1 1, r^) , (3.7) 



and the only allowed transition is fro m the state |1, — 1) 
to the state |1,0). Note that, in Eq. (3/7), we consider 
only the orbital part of the |a) states. 

Similarly we can introduce an analogous representation 
for the spin variable in the triplet spin states. This case 
is less straightforward, as we have to deal with the sum of 
two spins ^. This leads to a space whose dimensionality 
is four and we must project out the singlet subspace. 

The most direct way to introduce the spin represen- 
tation is to show how it is possible to write the corre- 
spondence between the matrix elements of the fermion 
and the spin operators, in analogy with Eq. ( |3.6[ ). Con- 
sidering all the possible spin transitions that leave us in 



the triplet subspace, \S ■ 
have: 



1, Sz), spanned by the \a)i, we 



Due to the properties of the Hubbard Hamiltonian, we 
have no transitions |1,1) ^^ |1,— 1). The labels m, m' 
and m" (=1,2,3) refer to orbital occupancy, according 
to Eq. ( |3.6| ). This implies that m ^ m! and ra ^ m" , 
but there is no constraint on m! and m" . The spin la- 
bel a (a) can be either t (i) or [ (f). Again, as in Eq. 
( |3.6[ ), such equivalence must be interpreted as equality 
of matrix elements between corresponding states in the 
I a) subspace. The only difference with the previous case 
is due to the factor -j= that appears each time a state 
involved in the transition has a component in the si ngle t 
state S = 0. For example, the second line of Eq. (3.S) 



(a|c+^c„xc+,^c™"t|a') = 



(3.9) 



and the only allowed transition is between triplet states 
with Sz = 0. Note that this time, due to the Pauli prin- 
ciple, we indicate e xplic itly t he o rbital labels tti, to', to", 
in contrast to Eqs. ( pl^ ) and (3/7), where we did not indi- 
cate the spin labels. Nonetheless, Eq. ( ^.g| ) expresses an 
equivalence between matrix elements in the spin space, 
only, the orbital degrees of freedom having already been 
taken care of. The global spin-orbit representation is ob- 
tained through the direct product of Eqs. (3.7) and (p!s|). 
For example, considering both the orbital and the spin 
degrees of freedom, we have the equality: 

(a|c+^cixc+j^C3||a') = (t = 1,Tz\ - -^t+Tz\t ^ 1,t^) 
X {S^ I, S zl'^il - S^,)\S ^ 1, S'J . (3.10) 

In this way it becomes possible to rewrite the effective 
Hamiltonian Hcs in terms of the spin and pseudospin 
operators. A straightforward but nevertheless tedious 
algebra leads to the expression for the spin-orbital iJoff 
which we report in Appendix C. Note that the main diffi- 
culty of this Hamiltonian consists in the great number of 
terms to deal with (3"* = 81) which is due to the fact that 
there is no conservation law for the pseudospin quantum 
number r^ in the hopping process. This prevents from 
having more symmetrical expressions for the orbital part 
of iJeff, in contrast with the spin part that retains the 
usual spherical symmetry. 



We can now make some general considerations about 
Hes- First it will be expedient to write it in the following 
form: 



H. 



■'« = -lu^ E,j [2 + 5^, • Sj] O 



(1) 



(2) 



4 (72+4J ^^i] 

J 

' 12 U2+2,J 



(3.11) 



WE.,[i-5.-5,]o^ 



(3) 



i(fc) 



where O^ are the orbital contributions to the energy 
corresponding to three different terms of the effective 
Hamiltonian H^s given in Appendix C. It is clear that 
the magnetic behavior of the system described by i^eff is 
strongly affected by the orbital degrees of freedom, that 
determine the sign and order of magnitude of the ex- 
change constants. Nonetheless, even without considering 
the explicit form of the O^^ , we can easily demonstrate 
t he fo llowing statement about i?eff starting from its form 
( |3.11 ): "It is impossible to get the ground state of each 
single bond in a configuration of total spin Si + Sj = 1" . 
The proof proceeds as follows. The spin scalar product 
can take the following values, depending on the kind of 
coupling between the two nearest neighbor spins: 



(3.12) 




^^7; ■ •^') — ^ ■ 



If we consider a single bond, ij, this implies the fol- 
lowing form for the Hamiltonian in the three different 
cases: 







where we defined: 



U.-J^'' 






(1) 



2 ttA 



fe-) 



and 



4C/2 



-47^^^ 



(2) 



4C/2 



-27^^^ 



(3) 



(3.13) 



(3.14) 



(3.15) 



ferromagnetic and the antiferromagnetic bonds, respec- 
tively. Note that the two Hamiltonians H^^ and H^ act 
as projectors on the subspaces of two spins coupled to 
Sij = 2 and Sij = 0. In the following we shall call them 
ferromagnetic and antiferromagnetic parts of -ffoff- 

Another important consequgnce is the impossibility to 
describe by iJeff the old CNRLl solutions, where the ver- 
tical molecule (see section V) was bound in a molecular 
spin Sij = 1, as a particular limit in the parameters 
space. The two models differ drastically from the begin- 
ning (the older being an S" = 1/2 atomic- limit descrip- 
tion), and there is no solution of the new model where the 
molecule is bound in a total spin Sij — 1 state, as long as 
the effective Hamiltonian HcS is a good representation of 
the problem. 

Finally we should remember to add to the second order 
HoS the first order contribution coming from the crystal 
field, whose form is easily seen to be 



rrc/ _ 
-"eff — 



EA. 



jz 



(3.16) 



0. 



where Ai = A_i = At and Aq 

Before ending this section, it might be instructive to 
present the result of a preliminary study that led us to 
the spin-orbital representation of Hes: the case of two 
electrons coupled to spin S = 1 and two orbitals per site. 
The usefulness of this solution is that one could derive 
the answer on physical grounds. In fact, due to the high 
J value, the crystal ground state is composed by a collec- 
tion of atomic states, each with one electron per orbital 
ferromagnetically coupled, in order to maximize Hund's 
exchange, as in the more general 3-orbital case. The 
main difference is that now there is no orbital degener- 
acy. Due to the "freezing" of configurations the only way 
to decrease the atomic energy with 2'"* order hoppings 
is to have an antiferromagnetic coupling between nearest 
neighbors with exchange constants given by an average 
over all possible atomic hoppings among the four orbitals 
involved. As expected, we obtained: 



Hr 



1 



iU2 + 3 J) 



J2{tW-S^■SJ], (3.17) 



It is easy to see that if H^{ij) < H^ff{ij), then the 
minimum of energy is achieved for the antiferromagnetic 
configuration Si + Sj = 0, while, if, on the contrary, 
H^ff{ij) < H^[ij) then the minimum of energy is for the 
ferromagnetic one Si + Sj = 2. Thus for any value of the 
parameters the minimum of the energy is never achieved 
in the configuration Si + Sj — 1, except for points in 
parameters space where H^ = H^g and all three con- 
figurations are degenerate. This gives us the important 
result that we can simply study the two spin configura- 
tions Si + Sj = 2 and Si + Sj — 0, corresponding to the the hexagonal ch axis (see Table II of next Section). 



where (if^) = Em,m'=i(^Ij" ? and U2 + 3J is the energy 
difference between the atomic ground state and the po- 
lar intermediate state. This derivation was also useful to 
find out the spin-1 representation in Eqs. ( t3.8[ ), (3^) in 
terms of fermion operators. 



The result in Eq. ( 3.17 ) confirms the expectation that 
with no orbital degeneracy it is impossible to break the 
trigonal symmetry in the corundum structure of V2O3, 
since (if,) is invariant under a rotation of 27r/3 around 



IV. THE CRYSTAL AND MAGNETIC 
STRUCTURE OF V2O3 



Before embarking in the study of V2O3, we shall try 
to establish the magnetic symmetry group of the AFI 
phase from the analysis of the present experimental evi- 
dence, since this will turn out to be important for decid- 
ing between different ground state solutions with the real 
spin structure obtained by the minimization procedure of 

The crystal structure of V2O3, the choice of the Wan- 
nier basis functions, their symmetry properties and the 
role of the oxyg^a, have been discussed in detail in the 
papers by CNR,Elfl to which we refer for what not ex- 
plicitly repeated in this paper. We adopt here the same 
conventions and definitions. For convenience we shall re- 
mind here some basic facts about this compound. 




monoclinic cell 
(primitive) 



Xj^sy and .^♦J*-' have "opposite" oxygen octahedra 

FIG. 1. Corundum structure together with the unit cells 
for the trigonal and monoclinic phase. Only Vanadium ions 
are shown; filled and empty circles correspond to ions with 
Oxygen octahedra differently oriented in the space. Arrows 
indicate the direction of spins. 

Figure shows the non primitive hexagonal unit cell 
of V2O3, together with the primitive trigonal one, in the 
corundum paramagnetic phase. Only Vanadium atoms 
are shown. There are two formula units per cell (four V 
atoms) and each metal ion is surrounded by a slightly 
distorted oxygen octahedron with the three-fold symme- 



try axis directed along the c/f-hexagonal vertical axis. 
The distortion corresponds to a compression of the oc- 
tahedron along cu- The octahedra around the V-ions 
represented by the filled circles in Fig. |l| are rotated by 
180° about the Cn axis with respect to those around the 
V-ions represented by the empty circles, this orientation 
varying from plane to plane. Only 2/3 of the octahe- 
dra centers are occupied by the metallic ion. The space 
group of the corundum phase is i?3c with the fplJowing 
generators (written in the conventional notationEj): 



{^,0} 
{4,0} 

{C2, o^i 



Identity 

Inversion around the midpoint of V4 
and V5 (point O in Fig. 1) 
Rotation of tt/S around the ch axis 
Rotation of tt around the &m axis 
followed by a translation of \am , 



where a^ and b^ are two of the basis vectors of the mon- 
oclinic unit cell shown in Fig. |l|. The origin O has been 
chosen as the inversion point between atom V4 and V5 
in the trigonal cell. The translation associated with the 
C2 rotation can be expressed also in terms of the trig- 
onal basis vectors defined in Fig. ^ as ^(a* + h + Ct). 
The corresponding crystal point group, obtained by set- 
ting all translations to zero, is Z^adBj Note that among 
the symmetry operations of R3c there is a glide plane, 
obtained through the combination of / and C2' 

{(Tf,, ^djn} = Refiection through a plane orthogonal 
to the bm axis followed by a translation 

of jflm . 

By lowering the temperature, the system makes a dis- 
ruptive first order transition to a monoclinic phase, with 
further distortion of the octahedra to accommodate a ro- 
tation by about 1.8° of the vertical Vanadium pairs (e.g., 
Vi and V4) in the a„i~Cm plane towards the adjacent oc- 
tahedral voidsJj As a consequence one bond in the basal 
plane becomes longer than the other two by about 0.1 
A, the trigonal symmetry is lost and the lattice space 
group lowers to I2/a with the same generators except 
for ((73,0). Its crystal point group is C2h- The mono- 
clinic cell is body-centered and, containing four formula 
units, it is not primitive, from the point of view of the 
bare crystal lattice. However concomitant to the struc- 
tural transition a magnetic order sets in, with ferromag- 
netic ttm " Cm planes stacked antiferromagnetically and 
with an AF wave vector given by |6rn-Q Because of this 
magnetic order, the monoclinic cell becomes primitive, 
due to the AF coupling of the magnetic moments on the 
V-ions connected by the body-centered translation. The 
magnetic moments of the V-ions, indicated by arrows in 
the figure, lie in the arTM-i— Cm plane, at an angle of 71° 
away from the ch axis.QQ Notice that the in-plane longer 
bond corresponds to the ferromagnetic coupling and is 
orthogonal to bm- 



Under these conditions one can easily check that the 
time-reversal operator T followed by the non primitive 



translation ^ [dm + i>„ 



is a symmetry operation of 



the magnetic structure, so that the magnetic point group 
is C2h ® T. Each operation should be followed by the ap- 
propriate translation as indicated here: 



l)E, I 
2) C2, db 

3) f , fi 

4) fC2, fab 



No translations 

2 y^m I ^7n) 



hi^r, 



-b„ 



c„ 



a„ 



Notice that the translation associated to C2 and ab 
has changed. In fact, since now the application of these 
two operations changes the direction of the magnetic 
moment ,e2I the total translation must be 






bm + ' 



— o y^rn ! ^7n) 



and the role of C2 and ab in the paramagnetic lattice is 
taken, in the AFI phase, by TC2 and Tdb- 

Under these operations the correspondence between 
the charge and magnetic states of the various metal sites 
with their oxygen environment is given in Table I, with 
reference to the numbering of Fig. m 

Table I. Correspondence table between magnetic sites in 
V2O3 



E : 


1 


2 


3 


4 


5 


6 


7 


8 


I : 


2 


1 


7 


5 


4 


8 


3 


6 


C2: 


8 


6 


5 


7 


3 


2 


4 


1 


&b ■■ 


6 


8 


4 


3 


7 


1 


5 


2 


T : 


7 


3 


2 


8 


6 


5 


1 


4 


TI: 


3 


7 


1 


6 


8 


4 


2 


5 


TC2 ■■ 


4 


5 


6 


1 


2 


3 


8 


7 


T&b : 


5 


4 


8 


2 


1 


7 


6 


3 



Now the recent observation of non reciprocal x-ray gy- 
rotropy by Goulon et allEl in the AFI phase of V2O3 
points to a reduction of magnetic symmetry. In this 
experiment a transverse x-ray linear dichroism at the 
Vanadium K-edge is observed and interpreted as due 
to a dipole-quadrupole interference effect. This signal 
changes sign according to whether the externally applied 
magnetic field is parallel or antiparallel to the direction 
of the incident x-ray beam along the ch axis. Therefore 
neither T nor / can be separately symmetry operations, 
but their product TI is. There are seven subgroups of 
four elements of the magnetic point group C2h®T. They 
are listed below: 

1. E,i,C2,ab 

2. E,i,f,fi 

3. E,i,fC2,fab 



E,T,C2,TC2 
E,T,ab,Tab 
E,C2,fi,f&b 
E,ab,fi,fC2 



It is immediately clear that only groups 6 and 7 are eligi- 
ble candidates, i.e., group C2h{Cs) and C2h{C2), respec- 
tively 2/m and 2/to in international notation. Both are 
magnetoelectric (ME); however the first one gives rise to 
an off-diagonal ME tensor whereas in tlic second one this 
tensor presents diagonal components Ed It is possible to 
discriminate between them by noting that the existence 
of an off-diaspnal ME tensor explains why Astrov and 
Al'shin failedEj to find a ME effect in V2O3, since their 
experiment was set up to look for diagonal components. 
This is a strong indication that C2h{Cs) is the correct 
magnetic group for V2O3. 

The origin of the reduction of magnetic symmetry from 
C2h 'S)T to C2h{Cs) can reasonably be ascribed to an or- 
bital ordering in the magnetic and charge density of V2O3 
due to electron correlatioas. However the ferro-orbital C 
phase found by Mila et altB does not provide the correct 
answer, since the corresponding magnetic group is easily 
seen to be C2/1 (^ T, due to the fact that all sites are oc- 
cupied by the same orbital. The same can be said for the 
other stable phases with the real spin structure found in 
this work. 

We speculate that the presence of an excited config- 
uration with the correct magnetic group C'2hiCs), very 
near the ground state and with the favorable coupling 
with the lattice, can provide the solution to this puzzle. 




FIG. 2. The neighbor structure of the vertical bonds. Re- 
ferring to Fig.[lJ, if the atom a coincides with Vi , then we have 
the correspondences 6 — > V4, c ^ V2, c' — *■ V5, d = e -^ V3, 
and d' = e' — > Ve. 

In order to proceed in the following sections with the 
minimization of HoS we need to have a reasonable guess 
at the various parameters appearing in it, namely the 
hopping integrals and the Coulomb and exchange atomic 
parameters. In Fig. we show half of the cluster of near- 
est neighbors to a given molecule (the other half can be 
deduced with the help of Fig. 1^) to illustrate the notation 
that will be used later on. 
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Following CNRul'B we present in Table II the transfer 
integrals evaluated by exploiting the symmetry proper- 
ties of the corundum structure without taking into ac- 
count the monoclinic distortion of the bonds, since again 
our aim is to show how electronic correlations can break 
the initial trigonal symmetry of the lattice. The devia- 
tions from this symmetry will be considered later on to 
illustrate if and how the monoclinic distortion can stabi- 
lize the orbitally ordered state with the correct magnetic 
spin structure. 

Table II. Transfer integrals along different bonds in the 
corundum phase. 



direction 


h 


h 


h 


Si 


til 


—a 


-\a+iP 


-\a+iP 


^^ 


^22 


/3 


-i« + i/3 


-ia+i/3 


Ai 


^33 


(J 


a 


a 


p 


il2 





^(« + /3) 


-%a + l3) 





tl3 





fr 


-^r 





hs 


— r 


kr 


kr 





t21 





^(« + /3) 


-^{a + P) 





hi 





fr 


~^r 





^32 


— r 


kr 


kr 






Table III. Transfer integrals (eV) from tight binding cal- 
culations used by CNRO and from LAPW-calculations by 
MattheissO 





Castellani et alB 


MattheissB 


M 


0.2 


0.2 


p 


-0.72 


-0.82 


—a 


-0.13 


-0.14 


P 


-0.04 


-0.05 


a 


0.05 


0.05 


— T 


-0.23 


-0.27 



In the following we shall assume for the Coulomb and 
exchange parameters j-the values suggested by Ezhov et 
ain and Mila et alB i.e., J ~ 0.7 ^ 1.0 eV, U2 ~ 
2.5 eV, and iar the hopping parameters those derived 
by MattheissH and shown in Table III. This set will 
be referred as the standard set. By fitting the LAPW 
band structure of V2O3 to a tight-band calculation, 
Mattheisa^ has provided the relevant Slater-Koster in- 
tegrals that have been used to calculate the appropriate 
hopping integrals. As on can see from Tf^lple III they are 
quite close to those estimated by CNRl 



V. THE ENERGETICS OF THE VERTICAL 
MOLECULE. 

As realized by CNRQ-eI and later by Mila et al.^ the 
formation of the vertical molecular bond {5^ in Fig. ||) is 



the key to the understanding of the physics of V2O3 in all 
three phases. This fact is indeed supportedrbjy the exper- 
imental evidence bothJiajm optical spectracj and inelas- 
tic neutron scattering.ErE3 However the solution proposed 
by CNR was appropriate to low values of J (< 0.2 eV) 
and values of the trigonal distortion which where sup^ 
posed to be quite small, as suggested by Rubinstein. E3 
The new solution proposed by Mila et al£^ reconciles 
the present evidence for an high value of J (> Q.l^s^fi 
and the consequent spin S' = 1 state of the V-iongl3'Efl 
with the existence of an orbitally degenerate molecular 
state, while being rather stable against a sizable value 
of the trigonal distortion. It is interesting to study how 
this can come about, since this investigation can provide 
a clue to the kind of variational wave function to be used 
in the minimization of ifcff, will delimit the regions of 
stability of the solution in the parameter space and in- 
dicate competing states that might be relevant for the 
phenomenon of the metal- insulator transition in V2O3. 



A. The approximate solution using H^s 



In considering the vertical pair it is convenient, as 
shown in Section V-B, to introduce the following molec- 
ular quantum numbers: the total spin 5*^ — Sa + Sb, 
its z-component S^^ = Saz + Shz and total z-componcnt 



of pseudospin t. 



M 



Tiyz- Along the vertical bond 



Si only ill — t22 = p and ^33 — p are different from 
zero (as seen from Table II) and their values are given in 
Table III. Specializing to this case the effective Hamilto- 
nian given in Appendix C we obtain for the ferromagnetic 
state {S^ = 2) the following expression: 



^Gff 



1 1 



3C/2- J 



[2 + Sa- Sb] ifl^f^ + p'fp 



P'PUp) , 
(5.1) 



and for the antiferromagnetic bond {Sm = 0): 
1 1 



■^cff 



2 C/2 + 4 J 



[l-Sa- Sb] 



1 1 



1 



X [^i {a^ + a'^)+ p ap + -ppa^p] 



1 — Sa ■ Sb] 



2U2 + 2J 



(5.2) 



X [p'^ia,, - a'p) + p^Gp - -ppa^p] 



[l-Sa- Sb] 



1 1 



6 C/2 + 2 J 



y^[2^? a"^ + 2p'a"p+iipa"pp] , 



where, respectively, we have defined: 
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Ju 



^2 _2 



'bz 



i - - + +_i 
2 ''a ''a ''"(, '''fc 2 



'a 'a 'b 'b 

Jp ~ '''az + '''bz 
Tp,p ~ '''a '^'azTbzT'i, + '''az^'a '''b '''bz 

+TaTaz'rbzT^ + TazTa '^b '^bz 



and 



a^ = 2 



21/ _l22n 

'^hz + -^^'''azTbz + '''az'''bz) 



O-'p = +4(^0 Ta ^6 Tfc 



'a 'a 'b 



-b+) 



A/jp 



Tq Tq2T^ T6^ + TazT^ Tb^T^ 



Ta '''azTfj T\,z 



1'azT^ TbzTij 



2 2 
az bz 



+ ''"62 '''azTbz 



2 2 
az bz 



+ T(^a ^a 7-b+T+ + T+T+T^ T^ ) 



'bz 



n2 Jl 



fJ-P 



Ta TazUzT^ + TazT^ T+Uz 
+TaTazTbzT^ + TazTa T^ Tbz 



We assume for the moment Af — 0. Based on Eq. 



(5.1) and with the definitions of Eq. ( |5.3| ) we can easily 
evaluate eigenvalues and eigenstates of H^^^. Neglecting 
the 5-fold spin degeneracy and taking into account only 
the orbital one, we find the following doubly degenerate 
ground state with rf^ = ±1: 



\r±)ab^ ^{\±l)a\Q)b+\±l)b\Q)a) 



(5.5) 



Equation (5.5) represents only the orbital part of the 
ground state. The whole state (e.g., with S^' = 2) can 
be pictured as: 



\l-)ab = -r. 



or as 




\l+)ab = -TS 



The corresponding ground state energy is: 

(M-P)' 



AE„ 



U2~J 



(5.6) 



With reference to the picture of the state, this energy 
lowering (with respect to the atomic limit 2(C/2 — J)) is 



made up of three terms: the virtual hopping back and 

2 
forth of an eg electron (— jj^_j ), the similar process for 

2 
the aig electron (— jj''_j ) and a sort of correlated hop- 

(5.3) ping in which an Cg electron jumps from atom a to atom 
b while simultaneously an aig electron jumps from atom 
b to atom a and vice- versa ( rr^^j , which is negative due 
to the opposite sign of p and /i, see Table III). This latter 
process is present only due to the "ent angled" orbital na- 
ture of the molecular state of Eq. ( |5.5| ) and is absent in its 

Hartree-Fock approximation, which provides a lowering 

2. 2 
of only AEhf = — ^ _ j . With the values given in Table 

III, the ratio between the interference and the HF term 
2iip/{^'^+p'^) is of the order of 50%. Therefore the molec- 
ular correlation energy AEm — AEhf = u^-J ^® much 
bigger than the in-plane exchange energy (« fj^n")? so 

(5.4) that in this case the best variational wave function for the 
entire crystal should be constructed in terms of molecular 
states. 

Another point worth mentioning here is the quality of 
the expansion around the atomic limit. The exact solu- 
tion of the 2x2 eigenvalue problem for the ferromagnetic 
vertical molecule is given by 



AE^ ^Era-2{U2- J) 



A{p-py 



iU2 Jf 



(5.7) 



We see from this expression that the expansion param- 



eter is 



2{p-p) 



,jj _js , which is of the order of one, using the 
standard values of Table III. This value is borderline for 
a good expansion; however, as often happens in pertur- 
bation theory, the second order term turns out to be a 
reasonable approximation to the exact result (1 -I- 0.89 
as compared to ^2.78 — 1.67, with a relative error of 
less than 13% in the worst of the cases). Moreover this 
problem is present only for the vertical pairs, since in the 
basal plane we are well within the values for a rapidly 
convergent expansion. Notice that, when comparing dif- 
ferent variational minimal solutions of Hcs in section VI, 
the error in the vertical pairs will cancel out and the re- 
sult will be of the same accuracy as the expansion for 
bonds in the basal plane. 

As long as \ p\ > \p\, the ferromagnetic state with 
5^^ = 2 in Eq. (5.5) is the ground state for the vertical 
pair. However it is easy to realize that, in the opposite 
case IpI < |/x|, the orbital part of the ground state changes 
to 



mab 



1 



V2 



(|l>a| 



(5.8) 



or, pictorially, including the spin (5* — 2) 



\l0)ab 
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with an energy lowering of AE'„^ ~ — tj^_ j ■ 

This state is not orbitally degenerate. However it is 
interesting to note that in this case the percentage of oc- 
cupation of the aig state is 50%, so that-this solution is 
excluded by the findings of Park et aL,Ej as well as on 
the basis of the theoretical estimates of Table III {p> fi). 

As long as ^ > 0.22 the ferromagnetic configuration 



in Eq. (5^) remains the ground state of the pair. By de- 
creasing J a transition to an antiferromagnetic (S"*^ = 0) 
ground state is expected, since this spin configuration will 
maximize the number of virtual hopping processes with- 
out loosing too much in Hund's energy. This is indeed 
what happens when -^ < 0.22. 

Even in this case we obtain a two-fold orbitally degen- 
erate ground state: 



1 



\4'^,\)abC^^{\-l)a\-l)b-\l)b\l)a) 



(5.9) 




and 



m 



1 

V2 



(|-1)„|1>5 + |-1).|1).), 



(5.10) 



mab^4^ 




\/2l 



For simplicity of presentation we have omitted to show 
the spin structure of this state, since this latter is given 
in th e foll owing section V-B (see states I71) and I72) of 
Eq. ( 5.17 )). The ground state energy is given by: 



AE. 



AF 



P 



2/i2 



U2 + 4 J U2 + 2 J 



(5.11) 



Note that the first state (Eq. ( |5.9| )), mixing the val- 
es r*^ = ±2 does not conserve the value of pseudospin 



r(4) 



z-component, due to the term Hq'^' in Eq. (2.2). 

The above level scheme is confirmed by the exact treat- 
ment of the vertical pair on the basis of the original 
Hubbard Hamiltonian which is reported in the follow- 
ing section V-B. There are slight discrepancies, however, 
due to the non optimal conditions for perturbation the- 
ory. For example, the transition value between the ferro- 
magnetic and antiferromagnetic configurations is found 
at -^ = 0.29. Moreover the degeneracy of the two anti- 
ferromagnetic states is removed, the second lying always 
lowest (from 1 to 5 meV, in the range of parameters of in- 
terest), due to the different mixing with states that have 



been projected out in the perturbation theory (those with 
Sa = Sb = and those with Sa = Sb = 1/2). In general, 
though, the exact energy level structure is reasonably 
close to the approximate one. 

By switching on the trigonal distortion At, the ferro- 
magnetic ground state energy ( |5.6| ) changes to 



U2 — J 



and the antiferromagnetic ( 5.11 ) becomes 



^^mtl - - 2 



2^2 



U2 -I- 4 J U2 + 2 J 



2At 



(5.12) 



(5.13) 



Because of this, the stability region of the ferromag- 
netic state in the parameter -d- initially increases with 
At, since its aig population is only 25%, compared to 
the 50% value of the antiferromagnetic states. However 
for At bigger than a critical value At there is an inversion 
of tendency and the stability region begins to decrease. 
This is due to the fact that the structure of the anti- 
ferromagnetic state changes abruptly from a situation 
in which the aig population is 50% to one in which is 
0%. Indeed for At > At the lowest energy configuration 
for the AF bond is reached when the two electrons on 
each site occupy both Cg orbitals and are coupled to spin 
5—1, with total spin S**^ = 0: in this case the ground 
state orbital configuration, for J not too close to zero 
(i.e., J > 0.2 eV, due to the constraint 2 J > At used in 
our perturbation theory) is given by: 



\ro^)ab ^ |0)a|0)6 = 




(5.14) 



The full spin structure of the state ( |5.14D will be given 
in the next subsection (see state I79)). The ground state 
AF energy, in this case, is: 



AEi'^, = -3m^ 



1 



1 



U2 + 4 J U2 + 2J 



(5.15) 



Fin al ly the esti mate o f At, obtained by comparing Eqs. 
( ^.12 ), (5.13) and (5.15), and using Mattheiss parameters 
of Table III, is At ~ 0.34 eV , in good agreement with the 
exact value calculated in section V-B (Af ~ 0.30 eV). 

We shall not dwell anymore on this subject since it will 
be studied more in depth in section V-B. 



B. The exact solution using the Hubbard 
Hamiltonian. 



In constructing the effective Hamiltonian (see section 
II) we have assumed that J is a high energy parame- 
ter and therefore have excluded from our zeroth order 
degenerate manifold singlet spin states, lying at higher 
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M 



energies by 2J or more. In order to assess the range va- 1) S 
lidity of i?eff as a function of J and the stabihty of the 
vertical molecule for J — > 0, we examine here the ground 
state configuration of the vertical molecule using Hub- 
bard Hamiltonian {Hq + H'). In this case the Hilbert 2) S 
space is made up of 495 atomic states with up to four 
electrons per site and the eigen-problem is rather com- 
plicated. However, due to the SU(2) invariance of the 3) S 
Hubbard Hamiltonian, states with different spin S**^ and 
different S^^ do not mix. Furthermore, since by assump- 
tion only diagonal hopping integrals are different from 4) S*^^ = 
zero, the kinetic and crystal field terms H' in Eq. ( |2.3| ) 
do not change the total molecular pseudospin r*^, while 
the atomic part Hq can only mix states with the same 5) 

r(4) 



parity of r^, due to the i?g term in Eq. (2.2). We 



can therefore divide the states of our Hilbert space into 6 
groups which are characterized by the following quantum 
numbers: 



^M 



7. 

M 



z 



M 



rM 
z 
M 



Z 

' z 



z 
M 



s 



rM 



--S^ = and T^ even (57 states): 

±4 (2 states); 

±2 (28 states); 

(27 states); 
: ^f = and T*^ odd (48 states): 

±3 (8 states); 

±1 (40 states); 
: 1, Sf fixed and r*^ even (49 states): 

±2 (22 states); 

(27 states); 

1, 5*=^ fixed and r^ odd (56 states): 

±3 (8 states); 

±1 (48 states); 
: 2, S^ fixed and r^ even (7 states): 

±2 (2 states); 

(5 states); 
: 2, S'f fixed and r*=^ odd (8 states): 

±1 (8 states). 
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FIG. 3. Energy eigenvalues of the Hubbard Hamiltonian for the vertical molecule as a function of J: solid line - S ' —2 
and T^^ odd; dot-dot-dashed line - 5*^^ — 1 and t^'' even; dashed line - S^^ — and t^^ even. Panels (a),(fe),(c) and (d) 
correspond to the trigonal distortion At = eV, At = 0.1 eV, At = 0.3 eV and At = 0.4 eV, respectively. 
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This classifications is valid for each S^ component, so 
that for S'*^ = 1 we get (49 + 56) x 3 = 315 states and 
for S"*^ = 2 we get (7 + 8) X 5 = 75 states. It is worth 



r(4) 



noticing that when the Hq'^' term in Eq. (2.2) is not 
effective, t^ is conserved and further reduction of the 
number of mixing states is possible. This is the case for 
all the states with 5^^ — 2, where the 15 states of fixed 
S^^ component can be grouped into the orthogonal sets: 
T^i = +2 (1 state), T^ = +1 (4 states), r*^ = (5 
states), t5 = -1 (4 states), r^ = -2 (1 state). 

We take the standard values of the parameters 
(Mattheiss set) and look at the eigenvalue structure as a 
function of J for different subgroups of the Hilbert space 
(see Fig. 0). For this choice of the parameters the three 
lowest energies are always in the groups of states with 



S'*' = and rf even, 5*^ = 1 and rf even or S^^ = 2 
and r*^ odd. They are presented in Fig. |[ 

The crystal field degenerate case (A^ = eV) is pre- 
sented in the Fig. ^(a). For J —>■ eV the ground state 
of the vertical pair belongs to the sector 5'*^ = with 
T^^ even. By increasing J three transitions occur. At 



very low J ~ 0.04 eV a first transition to a state with 
5*^ — 1 with T^ even takes place. This state^is remi- 
niscent of the ground state postulated by CNR0 for the 
vertical molecule (about 55% of its weight is composed 
by the old CNR state for J — 0.04) and, because of this, 
it does not belong to the Hilbert subspace upon which 
HcS can operate. In this region the 5*^ = with r*^ 
even state lies only 3 meV higher in energy. Then at 
J ~ 0.17 eV a second transition takes place again toward 
a state with S'^ = and t*^ even. To get an idea of its 
composition in terms of atomic states, we give its expres- 
sion at J — 0.73 eV, that is the upper boundary of the 
region of stability of the S^^ = state. Note, that even 
though the weight of the particular state depends on the 
value of J, the tendency of the weight distribution is the 
same in the whole region 0.17 eV< J < 0.73 eV. We get: 



||G^)j=o.73-0.67|7i)+0.67|72) 
+G.14|73) + 0.14|74) + 0.14|75) + O.U\-fe) , 
where 



(5.16) 



i7i; 



© 



© 



© 






1^^) = ^ 2 



l73) = ;^( 



l74)-;^l 



l75) = 4 



72 



l76) = :k 



V2l 








This state is essentially composed by 



G5)j=0.73^^(|7l) + l72)) 



(5.17) 



i.e., the t^ — combination of nonpolar atomic spin-1 
states coupled to 5**^ =0. Its orbital part is the same 
as the state |'0o)ab mentioned in Eq. (5.10), of which it 
is the complete spin orbital representation. Note that 
the same spin structure belongs also to the state |^['2|)a6 



given by Eq. (5.9), even if the orbital part is different. 

At still greater J, the value J ^ 0.73 marks the final 
transition to the doubly degenerate ferromagnetic state 
with S^'^ — 2 and t*^ = ±1, which is therefore stable for 
J/U2 > 0.29. We have two 2x2 eigenvalue equations for 
both T^ = ±1. Choosing t*^ = — 1 and solving for the 
ground state, we get: 



||G5)j>o.73 = A^(l77) + C|78) ) 



(5.18) 
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where C — -M^ 



U2-J 



(1^ V^ + '^wf)' 



2(/^-p) V V iU2-jy 

priate normalization factor and, for S**^ 



I77) = 4, 



N is an appro- 
= 2, 




l78> = :k 




V2I 




m\ ml\' 



Notice that the state in Eq. ( ^.5| ) is nothing else that 
the atomic limit of \GS)j>o.73 (actually, I77) = |7_)afc)- 
For the chosen values of the parameters and J = 0.73 eV 
the component I77) represents 83% of the total weight 
and increases with increasing J. 

To demonstrate the role of the trigonal distortion At 



on the stability region of the various ground states, we 
consider different values up-|to 0.4 eV, which is the value 
suggested by Ezhov et a/JlS As seen from Fig. ^(a) to 
Fig. |^(c), for values of At up to 0.3 eV the role of the 
trigonal splitting is essentially to decrease the value of J 
at which the 5*^ = ^ S^^ = 2 transition takes place, 
i.e., to increase the stability region of the ferromagnetic 
state. As already anticipated in the previous subsection, 
this fact can be easily explained by looking at the struc- 
ture of the S**^ — 2 state, essentially composed by I77) 
with 25 % of aio occupancy, and the S^^ = state, given 
in Eq. ( 5.17 ) and composed by I71) and I72), with 50% 
of aig occupancy. At At ~ 0.3 eV an abrupt transition 
in the composition of the S^' = ground state takes 
place such that the preferred orbital occupation change 
to the Eg states (no aig orbitals, see state I79)), while the 
S^'^ = 2 remains the same. As a consequence the sta- 
bility region of this latter starts decreasing, as shown in 
Figs. |^(c), (d). This transition in the composition of the 
S^'^ — state had to be expected, since for ■^, — — + cx3 
the aig occupation must go to zero. 
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FIG. 4. Energy level structure at J = 1.0 eV for molecular states with S" = 2, rj" odd; S'" = 1, r*' even and S" = 0, 
Tj even. Panels (a), (6), (c) and (d) correspond to the trigonal distortion At = eV, At = 0.1 eV, At = 0.3 eV and At — 0.4 
eV, respectively. Dot-dashed, dashed and solid lines indicate two, one or zero atomic aig orbital occupancy in the ground state. 
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Figure ^(a)-((i) illustrate how the low lying level struc- 
ture and the composition of the ground state of the ver- 
tical pair changes as a function of A^. We fix J = 1.0 eV, 
for reasons that will be apparent in the next section, al- 
though the same results are qualitatively valid for all the 
physical values of J. The main difference is that be- 
low the chosen value of 1.0 eV, the energy gap between 
the S^'^ — 2 ground state and the excited 5*^ = and 
5"*^ = 1 states reduces, while above J = 1.0 eV, it in- 
creases. At low values of At there are many low-lying 



l79> = 5^ 2 









excited levels for S^ = 0, t^ even states (dot-dashed 
lines of Fig.Q(a) and (6)). The orbital part of the two 
l owes t states in this situation is given by Eqs. (5.9) and 
( ^.10|) and their energy difference is about 6 m eV. N otice 
that the lowest state is exactly given by Eq. (5.17). For 
the chosen value of J = 1.0 eV, at At ^ 0.18 eV there is 
a drastic redistribution of the weight in the atomic con- 
figuration of the S**^ = molecule state, i.e., the non 
polar state is now given by 



m^ 





a^ 





which is nothing else th at th e complete spin-orbit rep- 
resentation of the state ( 5.14 ). By further increasing At 
(see Fig. p|(c) and {d)) this state becomes more and more 
favorable and finally at At ~ 0.47 becomes the overall 
ground state. At this value of At the weight of the non 
polar state [79) is more than 99%. Note that the value 
of the crystal field splitting where we have the transition 
from a S^'^ = 2 to a S**^ — overall ground state does 
depend on J as is clear from Figs. 0. 

For J — 1.0 eV the S**^ — 1, t^ even state is never 
the ground state of the vertical molecule. Nevertheless 
at At = 0.4 eV (see Fig. |((i)) it hes only about 80 meV 
above the ground state with S^ = 2 and its composition 
is made up for more than 99% of the following state (e.g., 
with Sf = 1): 



I710) 



As is clear from Fig. ^(d), its excitation energy de- 
creases with J. Note that I710) is made up of two atomic 
S = 1 states, so that it belongs to the subspace of Hcs- 

We shall make use of these findings later on in order 
to determine the various parameters in the AFI phase of 
V2O3. 



VI. THE MINIMIZATION PROCEDURE. 

In this section we look for all the possible orbital 
and magnetic ground-state configurations of the effective 
Hamiltonian by using a variational procedure. The trial 
wave function can be written in general as follows: 




where the state IV'^) refers to orbital occupancy and 10^) 
refers to spin occupancy on site n. In the following, 
we will use as a variational wave function jVPn) either 
an atomic state or a molecular one, with n labeling an 
atomic or a molecular site, respectively. 

Discarding for th e moment the single site crystal field 
part in Eq. (3.16) which will be easily dealt with, we 
observe that Hcs acts only onto two atomic sites at a 
time and factors into an orbital H°g and a spin H^g part. 
Therefore its average value over the above state takes the 
form: 



(*„|(*„|i/eff|*m)|*„) - 



(6.2) 



|vi/) = n„ |vi/„) = n„|^°)|0 



(6.1) 



Whereas orbital averaging in the first term will require 
some algebra, the second average in this equation, refer- 
ring to spin variables, is straightforward in a mean field 
treatment. For a ferromagnetic bond (5„ • Sm + '2)hf — 3 
and {Sn ■ Sm — 1)hf = 0, while for an antifcrromagnctic 

couphng, {Sn ■ S„i -I- 2)hf = 1 and {Sn ■ Sm - 1)hf = -2. 
As discussed in section V-A, the correlation energy 
of the ferromagnetic state of the vertical pair, defined 
as the difference between the exact ground state en- 
ergy and its Hartree-Fock approximation, is given by 
AEm — ^Ehf = 2/0/^/ (f/2 — J)- Therefore we can have 
two qualitatively different regimes of solutions: 
i) If this difference is much higher than the interaction 
in the basal plane, then the most appropriate variational 
wave function for the whole iJoff must be constructed 
in terms of molecular units, taking into account exactly 
the molecular bi ndi ng energy, with orbital wave functions 
\ij°) given by (^.5|). This means that the whole crys- 
tal consists of some ordered sequence of molecular units, 
whose internal energy is so high that it is energetically 
more favorable for the system not to break this structure. 
This seems to be the case for the values of parameters 
given in Table II. This state will be called the crystal 
"molecular" variational state. 
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ii) If instead it is the values of the exchange energy in 
the basal plane to be bigger than the correlation energy 
2/9yu/(C/2 — J), then the most natural variational wave 
function can be written in terms of single site atomic 
states, as will be shown in section VI-B. 



A. The case of the molecular variational state. 

As mentioned above, when the molecular correlation 
energy is much bigger than the in-plane e xch ange energy, 
then the best candidate for \ip'^) in Eq. (3.1) is given by 
a linear combi nat ion of the wave functions of the type 
shown in Eq. (5.5). 

Therefore for any molecular site n {ah, cd , dd! or ee' 
with reference to Fig. H) the orbital part of the trial 
molecular electronic wave function can be written as: 



|V'°) =cosV'„|V'°)„ +sinV'„|'0+)r 



(6.3) 



In order to construct the expectation value of Hcs over 
the crystal molecular wave function, we need to know the 
result of its application over states of the form (6^), for 
n and m molecular labels. Considering, for example, the 
two vertical pairs ah and cd of Fig. y], this state can be 
written as: 



\r.c')\rab) = 

[cos 1pab\i^-)ab +sinV'afc|V'+)afc) 
X ( COS V^cc' I V'- ) cc' + sin tAcc' I V'+ > cc' ) 



(6.4) 



The details of the calculations can be found in Ap- 
pendix D. Then the contribution coming from a ferro- 
magnetic bond along 5i is found to be: 



COS^ 4lab COS^ V'cc'Gi + (sin^ V'afc COS^ ilcc' 

+ cos^ il^ab sin^ ipcc')G2 + sin^ f/'afc sin^ ipcc'Gs 
+ sin 2'ipab sin 2ipcc'Gi , 

where 



402 



(6.5) 



Gi = 


1(2/32 + 2(72 + 4x2+4x2+402) 


G2 = 


1(2^2^2/32 + 2^2+4x2 + 2^2 


G3 = 


f(2a2^ 2^2+4x2 +4^2 + 402) 


G4 = 


fa/3 


u = - 


1 



(6.6) 



U2-J 



Note that the spin contribution has been already taken 
into account. We have also retained for future use the val- 
ues of the hopping integrals 9 — i^^ and x = t^^j which 
are zero in the corundum phascEl and can be different 
from zero in the monoclinic one. 

For the AF bond along the same direction we obtain: 



(*ab|(*cc'|^cfr(<5l)|*cc')|*a6) = 

C0S2 -0a6 C0S2 ^cc'Fi + (sin2 tp^b C0s2 tp^c' 

+ cos2 ipab sin2 ■ipcc')F2 + sin2 tpab sin2 Ipcc'Fs 
+ sin2ipabSin2iljcc'F4 , 

with the notations: 
^1 = ^ + f (4a2 + /32 + a2 + 2t2 + 4x2 + 402) 



(6.7) 



^(12a2 + 7/32 + 7^2 _^ ^4^2 ^ 20x2 + 2002) 

3 ' 4^ 



^2 = % + i(2a2 + 2/32 + a2 + 3t2 + Sx" + 302) 



+ ^(10a2 + 10/32 + 7cr2 ^ ;l7^2 ^ ;l9^2 ^ ;l702-) 

^3 = fii + |(a2 + 4/32 + a2 + 4t2 + 4x2 + 202) 

+ ^(7a2 ^ 12/32 + 7ct2 + 20r2 + 20x^ + 1402) 

(6.8) 



V - 
w 



U2+iJ 

1 

U2+2J 



Again, we already included the spin contribution. To 
evaluate the averages of iJeS along 82 and 5^, it is conve- 
nient to uae the invariance properties of the Hamiltonian 
(see CNRQ) under the trigonal D^^ symmetry. Perform- 
ing a C3 rotation around the vertical axis, the state along 
(5i is transformed in the corresponding one along 5^: 

G3|*dd')|*ab) = 

[cos(V'drf' - 27r/3)|*_),e' + sin(V'drf' - 27r/3)| *+),,-] 
X [cos(V'ah - 27r/3)|*_)ah + sin{i>ab - 27r/3)|*+)afc] . 

In the same way, applying C^ to the same state, we 
get the corresponding one along 82- Then the expecta- 
tion valu e of iJoff can be easily obtained directly from 
Eq. (U) and Eq. (|0|). 

By summing over the three nearest neighbors of the 
molecular site n in the horizontal plane (5i, 82, and (^3 
of Fig. 0), we can write the energy of the cluster in the 
compact form: 



(*|(*„|i?eff|*„)|*) = 
I]m=l,2.3[cOs2(V'n + 7m) C0s2(V'm + 7m) 

X {sGi + sFi) + (sin2(?/;„ + 7^) cos2(-0m + 7™) 
+ cos2(?/;„ +7m)sin2(V'm +7m))(sG2 +SF2) 
+ sin2(Vj„ + 7™) sin2(Vj„i + 7„i)(sG3 + sFs) 
+ sin 2(V'„ + 7m) sin 2('0m + 7™) (sG'4 + sF^)] , 



(6.9) 



adopting the following notations: s = 1 (s = 0) if the 
horizontal bond is ferromagnetic (F) and s = (s = 1) 
for the AF bond, 7™ — when m = 1 and jm = ±27r/3 
when m = 2, 3. 
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In the AFI phase the unit cell contains four of these 
clusters and eight V atoms. Therefore the energy per V 
atom, Ey , is given by the sum of these four energy con- 
tributions plus the four molecular energies AEmt of Eq 
(5.12), divided by eight. Referring to the atom number- 
ing of Fig. [^ we have the following expression for Ey'- 



Ey — ■§[{Ei2 + Ei3 + Ei^i) 
+ (£'37 + E27 + E27') 

-f (-E45 + ^46 + £^46') 

+ {Ee8+E58 + E58')] + h^AE^t 



(6.10) 
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where Enm is the appropriate energy of the horiz ontal 
bond nm, whether F or AF. Note that in Eq. ( |6.1C| ) each 
of the four terms {E12 + £'13 + £13')- (£37 + £27 + E27'), 
(£45 + EiG + £^46'), or (^68 + £58 + £58') is given by Eq. 
(p. 

The term \/S.Emt — — \ ^_j + i^t is constant with 
respect to the minimization angles because it is half the 
binding energy of the vertical molecule. For this reason 
in the following we shall consider E'y = Ey — ^AEmt, 
that represents the energy gain per V atom due to the 
intermolecular basal plane interactions with respect to 
this reference energy. 

We have then performed the numerical minimization 
of this expression with respect to all t he fo ur indepen- 
dent angular variational angles of Eq. ( |6.3[ ) in the unit 
cell, using the standard parameters as given in section. 
IV or reasonable variations around them. As in CNREI 
the minimizing angular values will provide the molecular 
orbital occupancy throughout the crystal. 

We have examined the following four magnetic phases: 

• AF phase - all three bonds 5i, 62 and S3 are anti- 
ferromagnetic; 

• RS phase - Si is ferromagnetic and §2 and S3 are 
antiferromagnetic; (this is the spin structure actu- 
ally observed in V2O3); 

• ARS phase - Si is antiferromagnetic and S2, S3 are 
ferromagnetic; 

• F phase - all three bonds Si, ^2 and S3 are ferro- 
magnetic. 



which correspond, respectively, to phases G, C, A, F in 
the work of Mila et al&i 

Figures g(a) and (b) show a plot of E'y as a function 
of J for all these magnetic phases. Notice that for fixed 
At, E'y depends only on the ratio J/U2 and scales like 
T^/[/2, if T is the largest hopping integral in the basal 
plane. 
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FIG. 5. The energy gain E'y per V atom as a function 
of J for different spin configurations (AF, RS, ARS and F 
type), represented by dashed, solid, dot-dashed and dotted 
lines, respectively. In panel (a) the hopping parameters are 
the standard ones (Mattheiss set), while in panel (b) the set 
considered by Mila et al. is used (r = 0.27 and all others 
hopping integrals equal to zero.) 

One general feature that is apparent from these fig- 
ures and the next Fig. is that the stability region for 
the RS phase is very much reduced in the two param- 
eter space of the hopping integrals ancLthe ratio J/C/2, 
in contrast with the spin S = 1/2 case.Q The AF and F 
phases occupy nearly all phase space in such a way that 
J/C/2 ~ 0.4 almost marks the transition from a stable 
antiferromagnetic in-plane spin structure to a ferromag- 
netic one (see, for example. Fig. |^(a)). This behavior 
depends on the fact that for low J values, the system 
wants to maximize the number of electron jumps (which 
occur more easily in an AF structure) to the detriment 
of on-site Hund's energy gain, whereas for high J val- 
ues this last mechanism is prevalent. In between, in a 
small range of J /U2 slightly depending on the values of 
the hopping integrals, typically 0.30 < JIU2 < 0.45, find 
their place the RS and the ARS phases, each occupying 
about half of the interval. Even in the most favorable 
case (see Fig. ^(6)), their stabihzation energy, due to 
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the competing presence of the AF and F phases, is very 
small, of the order of 2 meV, to be compared with a 
transition temperature corresponding to 15 meV. Even 
though the stabilization energy scales like r^/C/2, there 
is not enough room to improve substantially the situation 
by a reasonable variation of the parameters. 

More specifically, we see from Fig. ||(a) that assum- 
ing Mattheiss' parameters for the hopping integrals we 
achieve a stable solution for the RS structure only for a 
very small window of the ratio J/U2 around 0.4. The 
most favorable situation for thisJatter is obtained with 
the choice made by Mila et aZ.jl3 by putting a — (3 — 
a — and r ^ (Fig. ||(6)), but even in this case, as 
already stated, the stabilization energy is of the order of 
K, 2 meV. We have also tried to investigate the role of the 
monoclinic distortion in stabilizing the RS structure. To 
introduce it, we have assumed that after the setting in 
of the broken symmetry phase, the hopping parameters 
X = i^^ and 9 — t^^ take a value different from zero along 
the bond Ji, whose length increases by 0.1 A, and remain 
substantially zero along the other two directions, where 
the bond distance is unchanged after the transition. The 
result is essentially negative as there is a little but not 
significant improvement. 

Turning now to the orbital structure. Fig. P(a) gives 
the mini mizing values of the orbital mixing angles %pn of 
Eq. (6.3) in the basal plane at all molecular sites for the 
RS configurations found in Fig. |^. The value ijjn ~ 
for all sites means a uniform occupation throughout the 
crystal of the molecular statG-41'0-)i i-^-, the ferro-orbital 
solution found by Mila et a/.Ej In agreement with them, 
also for the ARS configuration (their A phase) we find 
a uniform solution with a minimizing angle ipn — 7r/2, 
i.e., a uniform occupation throughout the crystal of the 
molecular state ||V'+)- 

Moreover for the AF and F phases we find a contin- 
uum of orbital degeneracies of antiferro-orbital type, in 
the sense that all the orbital configurations with any mix- 
ing angle ipa on the central molecule and a mixing angle 
of V'a + 7''/2 on the three in-plane neighboring molecules, 
have the same energy (see Fig. H(6) for the particular 
case ^a = 0). This feature can also be deduced analyti- 
cally from Eq. (6.9). 

As anticipated in Section IV, the magnetic group of 
the ferro-orbital solution is not in keeping with the ex- 
perimental findings of Goulon et aLtil We have therefore 
analyzed the orbital order of the excited configurations 
within a range of ~ 4 meV from the ground state. Refer- 
ring to Fig. ^(&), for J = 0.85 meV, the ferro-orbital 
RS(FO) phase is at E'y = -0.0744 meV, all the de- 
generate AF(AO) phases with antiferro-orbital ordering 
at E'y — —0.0729 meV and a phase that will be called 
RS(ME), for reasons that will shortly become apparent 
(ME stands for magneto-electric) , with the orbital order- 
ing depicted in Fig. ||(c) lies at E'y = -0.0707. 

We have also explored the consequences of varying the 
ratio a/r from the zero value assumed in Fig, uib) to 
about one. From the picture following Eq. (5.5) it is in 



fact evident that an a = in value of the same order of 
T = ^23 would favor an antiferro-orbital coupling along 
the 5i bond in the RS phase. 
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FIG. 6. The angles indicate the orbital occupancy of the 
vertical molecule at a given site, (a) Orbital ordering in the 
RS(FO) phase. (6) Orbital ordering in the RS(AO) phase, (c) 
Orbital ordering in the RS(ME) excited phase. 
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J/Uj 

FIG. 7. Phase diagram in the {a/r, J/U2) parameter 
plane. Here AF, F, ARS denote the corresponding type of 
magnetic order; RS and RS' denote the RS(FO) and RS(AO) 
phase with orbital structures shown in Fig. M. The solid 
lines indicate the phase boundaries. We put for simplicity 
/3 = o- = 0. 

To this purpose we have drawn the phase diagram for 
the various magnetic configurations in the plane a/r ver- 
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sus J/U2- The relevant result is shown in Fig. 0. As 
expected another RS' phase with the same spin config- 
uration as the RS phase appears, centered around the 
values a/r — 0.7 and J/U2 — 0.40, with the molecular 
orbital ordering depicted in Fig. 0(6) (in-plane antiferro- 
orbital (AO) ordering). Again an energy analysis of the 
excited phases for J = 1.0 eV and a — —0.19 eV leads 
to the following sequency: 

E'y[RS{AO)] = -0.0909 eV 
E'y[AF{AO)] = -0.0908 eV 
E'y[RS{ME)] = -0.0905 eV 

where in brackets we have indicated the spin and or- 
bital configurations. As seen, the energy spreading is now 
much more reduced, only 0.4 meV separating the orbital 
ME phase from the ground state. 

In order to establish the magnetic group for the RS 
phases, of interest here, we observe that, because of the 
entangled nature of the molecular state, the average or- 
bital type of the two atoms constituting the vertical pair 
is the same. With respect to the I?3d corundum sym- 
metry point group, the state |0) transforms according to 
the totally symmetric representation {Aig in Schoenflies 
notation) while | — 1) and |1) are partner functions of the 
bidimcnsional Eg representation and transform, respec- 
tively, like the basis e^^^ and ef^ in CNrII 

We can then use Table I to see which symmetry oper- 
ations conserve the colored magnetic structure, obtained 
by adding to the lattice sites not only a spin label but 
also a color label given by the type of orbital occupa- 
tion at that site. We find that the magnetic group for 
the RS(AO) phase is C2 «) T (group n. 4 in Section IV), 
whereas that for the RS(ME) phase is C2h{Cs) (group n. 
7). Therefore the ouly phase compatible with the find- 
ings of Goulon et aLcil is this latter. As it will be argued 
in Section VII, it might be possible that the combined 
effect of the symmetry breaking of the spin and orbital 
degrees of freedom lead to a favorable coupling to the lat- 
tice in such a way as to stabilize the RS(ME) phase with 
respect to the competing configurations. In this case the 
role of the monoclinic distortion would be essential to 
achieve the ground state with the correct symmetry. 



B. The case of the atomic variational state. 



Even though the values of the hopping parameters as 
shown in Table III seem to favor what we called the 
molecular regime, nonetheless we think it could be useful 
to analyze, for the sake of completeness, also the other 
regime. In this case the orbital on-site part of the varia- 
tional wave function is atomic-like and should be written 
as 

II i:°) = cos6'i|0)i +sin6'i(cos-(/'i|l)i -l-sinV'il - l)i) , 



allowing all the three states |0)i, |l)i and | — l)i to be 
present without any a priori restriction on their rela- 
tive weight, contrary to the molecular case. The relative 
weight of the three states is then determined through the 
minimization procedure with respect to the variational 
parameters 9i and ipi. The only restriction we impose is 
that the solutions must fill the whole crystal, with a pe- 
riodicity not less than that of the monoclinic cell. This is 
indeed a quite reasonable request, as the solutions with 
periodicity of more than the unit cell should describe ex- 
cited states. As the unit cell of V2O3 is formed by 8 V 
atoms, there are in principle 16 minimization angles. In 
order to simplify the problem, we use the symmetry re- 
lations between the variational angles dictated by all the 
possible magnetic space groups for V2O3 described in sec- 
tion IV. Indeed for each group, the states of the V atoms 
inside the cell are not independent, but are related by the 
symmetry operations, thus providing a reduction of the 
number of parameters. By taking the absolute minimum 
we shall determine the orbital and spin nature of the 
ground state, together with the corresponding magnetic 
group. In this way we exclude solutions not invariant 
with respect to the chosen groups, but we note that all 
the interesting subgroups for V2O3 have been taken into 
account. 

Table IV. Number of independent variational angles in 
the unit cell according to the various possible magnetic 
groups. 
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In Table IV we list, for each magnetic group, the 
number of independent angles associated with the cor- 
responding atomic sites. This number is obviously given 
by 16 divided the order of the group. 

Given the full expression for Hcs reported in Appendix 
C, we can then evaluate the matrix elements for H^^f and 
H^ (as defined by Eqs. (3.14) and ( 3.15| )) along vertical 
{64) and horizontal {51,62,63) bonds. The spin averages 
are again calculated and included in the formulas as in 
the previous subsection. 



In the case of H^g the matrix element for the ^4 bond 



is: 
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(cos^ 9a sin^ 9b + sin^ 9a cos^ 9b)Gi 



>■,!, o^n O-i/., r7„ 



(6.11) 



where 



sin^ 0a sin^ ^b sin 2%pa sin ^tptG^ 

sin 20Q sin 26'h (cos V'a cos ipb + sin ■0a sin i}jb)G4 , 



Gs = 2m^2 
G3 = -u^2 
G4 — —u^p 

U^-JJ^ ■ 

For H^ it can be written as 

.9/1 .9/1 7^ 



(6.12) 



Vb = 



COS^ 0a COS^ 9bFl 

.,2 a „;„2fl, L sin^ 61a cos2 6*6) A 



+ (cos' 



I sin^ t/6 



- sin^ 6a sin^ ^b-Fs 

- sin^ 0a sin^ 0f,(cos^ Va cos^ ^b + sin^ V-'a sin^ ipb)F4 

l^ 0a Sin^ 0b(cOS^ 0a sin^ 06 + Sin^ Va COS^ 0'6)-F5 
0„ sin 9i, sin 27/;„ sin 2?/;j. Fk 



(6.13) 



+ sin 
+ sin 
+ sin 20a sin 



I i^a •3111 lyfti^i^vjo i/^a '3111 v-^0 I 

I 0a sin 0b sin 20^0 sin 2'0bi^6 
"" sin 20b (cos V'a cos '0b + sin^Aa sin-0b)i^7 



where 



Fi =2(w + w)Ai2 
F2 = (i; + ^)m2 

F3 == (u + U))p^ 



¥p^ 



iw ,,2 



F5 = ^^i^ 

^6 = (f-f)/i' 

^7 = (§ - f W 



(6.14) 



1) ^^ ^^ 

(72 +4 J 



"^ ~ U2+2J ■ 

For the horizontal bond Si , the -ff^ average value takes 
e form: 



{^,\{^a\H!'^{Sl)\^a)\^c) 

= [sin^ 0a cos^ 0a(cos^ 0c + sin^ 0c sin^ 0c) 
+ sin^ 0c cos^ Vc(cos^ 9a + sin^ 0a sin^ Va)]^^! 

■^^ .0.. f-i_ .o„ o> 



+ [sin^ 0a sin^ 0'a(cos^ 0c + sin^ 0c cos^ Vc) 

+ sin^ 0c sin^ 0'c(cos^ 0a + sin^ 9a cos^ l/'a)]<52 

+ [2 COS^ 9a cos^ 0c + sin^ 0q COS^ 0a cos^ 0c 

+ cos^ 0a sin^ 0c cos^ Vc + sin^ 0a sin^ 0c(sin^ Va 

+ sin ipc) — sin 20a sin 20c sin ipa sin 0^2] G3 

+ sin^ 0a sin^ 9c sin 2-0a sin 20cG4 

+ sin 20a sin 20c cos i/ia cos 0cG5 

+ (sin 20a sin^ 0c cos ■0a sin 2-0c 

+ sin 20c sin^ 9a cos 0c sin 2-0a)G6 

+ sin 20a sin 20c sin ■0a sin tpcGj 

+ (cos^ 9a sin^ 0c + cos^ 0c sin^ 0a)G8 

+(cos^ 9a sin 20c sin tpc 

+ cos'^ 0c sin 20a sin 0a ) Gg 

+ [sin 20a sin i/;a(cos^ 0c + sin'^ 0c cos 20c) 

+ sin 20c sin i/;c(cos^ 9a + sin^ 0q cos 2^a)]Gio , 



(6.15) 



where 



Gi — uc? 
G2 - M/32 
Ga - ur^ 
G4 = ua(i 
G5 = uacr 
Ge — —uar 
Gr = —uPa 
Gs = ua"^ 
Gg = —uar 
Gio == uPt . 



For the same bond in the case of H^ we obtain: 



(6.16) 



(^c|(*a|Hc1f('5_l)l*a)|*c) = 
COS^ 9a COS^ 0c Fi 



+ (cos^ 9a sin^ 0c sin^ 0c + sin^ 0a cos^ 0c sin^ ■0a)F2 
+ (cos2 0a sin^ 0c + sin^ 0a cos^ 0c)F3 
+ sin^ 0a sin^ 0c(cos^ 0a sin^ Vc + sin^ ipa cos^ 0c)F4 
+ sin^ 0a sin^ 0c F5 

+ Sin^ 0a Sin^ 0c(cOS^ 0a + COS^ 0c)F6 

+ sin'^ 0a sin^ 0c sin 2^0a sin 20c F7 
+ sin 20a sin 20c sin 0a sin ■0cF8 
+ (sin^ 0a sin 20c sin 20a cos 0c 
+ sin'^ 0c sin 20a sin 2^pc cos ipa)Fg 
+ sin 20a sin 20c cos ipa cos ipcFio 

|-[sin20aSin?/'a(l + sin^ 0cSin^^0c) 

'-sin20cSinV'c(l-^''^'^^^ "^"^ ''' ^^ 



(6.17) 



I Loiii ^i/a -^iii v^a VI- I '^111 i/c '3111 Ycj 

|-sin20cSin?/'c(l +sin^0a sin^ 0a ) ] A 1 
I- [sin 20a sin ■0a ( 2 — sin^ 0c ) 
■' sin 20c sin^0c(2 - sin^ 0a)]F'i; 



^12 , 



with the definitions 

Fi = via" + /32) + f (3(a2 + p"") + Sr^) 
F2 = w(-a2 ^ ^2) _^ 1(^2 _ ^2 _ 2t2) 

F3 = v(a2 + r2) + H(2a2 + 3^2 ^ 2a2 + Sr^) 

F4 = v(a2+/32) + |(a2+/32) 



^5 = 


■(;a2 + f(3a2 


+ a2 + 5r2) 


^6 = 


^2 + f (-3a 


2 + 3/?2 + 2t2 


F7 = 


(-! + !)«/? 




F8 = 


f(/3a + T2)- 


f(/3a + r2) 


h = 


-(! - f W 




F10-- 


= -(f-fW 




Fn-- 


= -f/3r 




F12-- 


.-far. 





(6.18) 



The matrix elements along 82 and ^3 are easily derived 
from the expressions (3.15)-( |6.18 ) using the symmetry 
properties of HcS under the trigonal G3 symmetry, as 
done previously in the molecular case. 

In this way the average of the Hamiltonian on a fer- 
romagnetic bond is given by the sole H^^ contribution, 
while for an antiferromagnetic bond we have to add \F[^g 
to the contribution of H^^. 
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The term due to the trigonal field splitting (3.16) can 
be also taken into account. Its energy contribution per 
V-atom is given by: -^ J2i=i ^^"^^ ^i- 

Given all these ingredients, we can now evaluate the 
ground state energy for all 8 groups: the "true" ground 
state energy per V atom, Eat^ is then obtained as the 
absolute minimum among the 8 minima, determining in 
this way also the magnetic group for V2O3. The results 
are presented for two choices of hopping parameters cor- 
responding to those of Fig. §(a), ib) in the molecular 
regime and p — —0.82 and fi — 0.2 eV (Mattheiss set, 
see Table III) . We consider for the moment only the case 
At = 0. As apparent from Fig. H(a), (6), we have in 
the atomic regime one more curve (labeled VAF) giving 
the ground state energy of the system in a magnetic con- 
figuration in which the two atoms along 5^ are coupled 
antiferromagnetically, independently of all the other spin 
couplings along ^i, 82 and ^3. Of course, such ground 
state configuration was absent in the molecular regime, 
where we started from a S**^ — 2 state. It relates to the 
S^^ = molecular solution in Fig. y. 

Notice that a direct comparison between Fig. g and 
Fig. H may be misleading, since in the atomic regime the 
vertical bond (54 is included in the ground state energy 
Eat, while in the molecular regime we have subtracted 
from Ey the binding energy of the molecule: "^"' . This 
means that the energy ^f^ (~ -0.33 eV , for J ~ 1.0)) 
must be added to E'y of Fig. ||, thus restoring the correct 
numerical correspondence between the two cases. 

In particular, for values of J such that the VAF is not 
the stable phase, we can write the atomic ground state 
energy per V atom as 



E„, 



E'at 



2 C/2- J 



(6.19) 



where we separated the in-plane contribution E'^f from 
the energy gain along 84, i.e. 



2 U2-J ■ 



In this way, dropping all the common Hartree-Fock 
terms along the vertical bond ^4, we can write the two 
inequalities: 



Kt < E'y 

E't > E[, 



U2-J 



(6.20) 



The reaso n for these inequalities is the following. The 
first of Eqs. ( ^.2C| ) says that the in-plane energy gain with 
the atomic variational wave functions is always lower 
than the corresponding molecular one. This is to be ex- 
pected, since the variational space in the molecular case is 
reduced with respect to the atomic one, where the states 
are not constraint to satisfy the form of Eq. (5.5). 

The second of Eqs. (6.20) states, instead, that the 
molecular ground state lies lower than the atomic one, 
because of the correlation energy. Note that, while the 
first Equation is always valid, the validity of the second 



is limited to sufficiently high values of p and p (for ex- 
ample, the standard set) and its breakdown marks the 
transition point between the molecular and the atomic 
regime. 
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FIG. 8. The ground state energy Eat as a function of J 
for different magnetic configurations. The hopping parame- 
ters are taken as follows: (a) p — —0.82, p — 0.2, a — 0.14, 
13 = -0.05, a = 0.05, r = 0.27; (b) p = -0.82, p = 0.2, 
a = 0, /3 = 0, cr = 0, r = 0.27; (c) p = -0.82, /i = 0, 
a = 0.14, /? = -0.05, a = 0.05, r = 0.27. 
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The analysis of the two Figs. H(a), {b) shows that the 
RS phase is reahzed in both cases (though in Fig. ||(6) it 
is very smah). Indeed in the case of the standard set of 
parameters, the atomic solution is even more stable than 
the corresponding molecular one, the first excited state 
lying ~ 3 meV above. Nonetheless it cannot be consid- 
ered a good grou nd st ate solution for the reasons stated 
above (see Eqs. ( 3.201 )). Moreover, the magnetic space 
group of the solution is not ME, being group 2. of Table 
IV. The orbital pattern is made of planes of V-atoms, 
orthogonal to the c/^-axis, in |0) state, alternating with 
planes of V atoms in | — 1) state. 

Figure ^{c) presents the case p = —0.82 eV and fi — 
eV, in order to satisfy the criterion for using an atomic 
variational function, i.e., the correlation energy less than 
the in- plane exchange energy. In fact, in this c ase t he 
condition /i = implies that the second of Eqs. ( |6.2C| ) is 
not satisfied, and then Eat < Ey- Thus the RS solution 
obtained in this case is the overall ground state of the 
system for this value of the parameters. Unfortunately, 
this solution has the usual drawbacks already analyzed 
(stability energy of only 3 meV, space magnetic group 
2. of Table IV, not ME) and moreover the ME solution 
lies very far from the ground state (~ 100 meV), thus 
confirming the fact that the choice /i = is not suitable 
to describe V2O3. 

The effect of the crystal field in this situation is to fa- 
vor the occupancy of the |0) states on all the atoms, as 
expected. In the case of high trigonal distortion At = 0.4 
eV, we obtain (not shown) a direct transition from the 
VAF to the F phase, with no stable RS solution. In VAF 
phase, all the V atoms are in the |0) configuration, while 
in the F phase, the percentage of |0) states lowers to 
85%, the remaining 15% being essentially |1), to allow 
the hopping process. 

Note that even in the case of such a high value of the 
trigonal distortion (At = 0.4 eV), we checked that the 
molecular phase is by far the more stable state (for ex- 
ample, for the standard set of parameters and J = 1.0 
eV, the molecular ground state lies 60 meV lower than 
the atomic one). 



VII. DISCUSSION AND CONCLUSIONS. 



which the pccupancy of the aig orbital is 25%. However 
Park et alS3 report an occupancy of 17% for this orbital 
in the AFI phase. This is an indication that other states 
with only Cg occupancy are mixed in the ground state 
of th is phase. In reality, even though the state in Eq. 
(|6.4[ ) is the main component of the wave function, in do- 
ing the variational procedure we have neglected excited 
states of the vertical molecule lying within the range of 
the exchange energy in the basal plane. Now, since for 
all the RS phases the ratio J/U2 is around 0.35 -;- 0.40, 
our study of the vertical molecule in section V-B (see Fig. 
(Q)) shows that we should take At around 0.4 eV if we 
want to stabilize the ferromagnetic coupling and find ex- 
cited states with only Cg occupancy close enough to the 
ground state. Fig. ^d) illustrates the level structure in 
the case J = 1.0 eV, from which one can infer that two 
more states with S^^ — 1 and S^^ — can be mixed with 
the S^^ ~ 2 wave function. Notice from Fig. a{d) that for 
J — 0.85 eV these states are even closer to the 5*^ = 2 
ground state. A more general calculation will be done in 
the future. We can however estimate in an approximate 
way the amount of mixing by writing the variational state 
I*) = a\S^'^ = 2) + (3\S^^ = 1)+ -f\S^^ = 0) and impos- 
ing the condition that the average value of the molecular 
spin be 1.7, as derived from neutrons and non magnetic 
resonant scattering data (see Introduction). Therefore 

2a^ + lf3^ + 0j^ ^ 1.7 , 

which, together with the normalization condition 

a'^ + P^ + "f'^ — 1, gives 

^2 = 0.85-/3^2 ; 7^ = 0.15 - /3V2 • 

In the situation shown in Fig. 0(c), even the small spin- 
orbit interaction of the AFI phase of V2O3 (~ 25 meV), 
which couples only to S^ = 1,2 states and not with 
S"*^ — state, can make the ground state energies for 
S^'^ = 1 and S^ = comparable. This means that we 
can reasonably take (3'^ =7^, thus obtaining a^ — 0.8 
and 0^ ~ ^"^ ~ 0.1. Therefore the aig orbital occupancy 
reduces t0|_Q.8 x 25% — 20%, in good agreement with 
Park et a/o findings. This means that also Park's find- 
ings point to a value of the trigonal distortion around 0.4 
eV. 



In this concluding section we shall review the implica- 
tions of the spin 5 = 1 model for V2O3 described above 
in relation to the present experimental evidence. The 
comparison will allow us to focus on merits and draw- 
backs of the various solutions obtained in the previous 
section. 



A. Occupancy of a\g orbital and spin. 

In the molecular regime, the variational wave function 
was assumed to be of the form given in Eq. (p^), in 



B. Trigonal distortion. 



^rI 



In CNRH the trigonal distortion At in the corundum 
phase (and therefore also in the monoclinic phase, since 
the cage of the oxygens rerriains almost unaltered at the 
metal to insulator transition^) was assumed to be negligi- 
ble, in keeping with the suggestion by Rubistein.Eil This 
latter was based on the following experimental evidence 
in the paramagnetic phases of V2O3: 

• There is no detectable anisotropy in the Knight 
shift of the NMR spectrum; 
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• There is no observable quadrupole splitting in the 
NMR spectrum. 

• Zero anisotropy was observed in single-crystal sus- 
ceptibility measurements for various crystalline ori- 
entations with respect to the magnetic field. A 
large value for At would result in highly anisotropic 
Van Vleck susceptibility. 

However the theoretical estimate by Ezhov et al&B and 
the previous discussion on the aig occupancy seem to 
point out to a value of At around 0.4 eV. The problem 
is how to reconcile this high value with the above ex- 
perimental findings. The concept of an entangled state 
for the vertical pair described by the wave function in 
Eq. (^) solves this apparent inconsistency in an ele- 
gant way. Indeed, the short range magnetic, fluctuations 
observed by neutron scattering experimentsO in the para- 
magnetic phases imply a statistical occupancy of the two 
orbitally degenerate npJecular states of Eq. (5.5) and x- 
ray absorption spectraEJ show evidence of an average 25% 
occupancy of the aig and 75% of Cg orbitals. This is very 
near to the totally symmetric occupation of the three 
orbitals (1/3 each), restoring in such a way the cubic 
symmetry. Notice that in a static picture the conclusion 
of a negligible trigonal distortion would be unescapable. 



C. Monoclinic distortion. 

In their work Mila et altB point out that the ferro- 
orbital molecular ordering found for the C phase (our 
RS(FO)) causes an effective uniaxial stress on the lattice 
degrees of freedom, leading to a uniform rotation of the 
vertical V-V pairs, in agreement with the raonoclinic dis- 
tortion proposed by Dernier and Marezio.H At first sight 
this explanation seems convincing. However on second 
thought it is not clear why atoms of two adjacent ver- 
tical pairs belonging to the same basal plane (e.g., Vi 
and V2, or V4 and V5 with reference to Fig. m) should 
go on opposite sides of the line joining them (x-axis in 
the figure). Since the two V atoms with their local oxy- 
gen environment are on average in the same electronic 
state, the initial displacive force should have the same 
sign along the x-axis, by an application of the Feynman 
theorem to an Hamiltonian depending on parameters (in 
this case the position coordinates of the V atoms). While 
a posteriori we see that the lattice energy is lowered if 
the two atoms move toward the oxygen voids, the initial 
drive should go on the correct direction. This is indeed 
what happens if one assume the orbital ordering found in 
the RS(ME) phase. The two atoms in the basal plane are 
in different electronic states, one odd with respect to the 
operation x -^ — x (state |V'-)) and the other even (state 
|V'+)), and this situation is reproduced for each pair of 
atoms in the basal planes throughout the crystal. 

We observe that in the RS(ME) phase the stress aris- 
ing from the magnetostrictive forces due to the broken 



symmetry of the spin degrees of freedom and the one 
originating from the orbital coupling to the lattice, pre- 
senting the same pattern of broken symmetry, act along 
the same axis, in contrast to the other two phases. This 
fact might give an advantage to the RS(ME) spin-orbital 
configuration in its coupling to the lattice so as to become 
stabilized with respect to them. In such an instance the 
monoclinic distortion would be essential to achieve the 
ground state with the correct symmetry indicated by the 
Goulon experiment. Clearly more investigation is needed 
on this point. 



D. Anomalous diffraction. 

In their paper Mila et alS3 claim that the ferro-orbital 
ordering they projjose, while different from that invoked 
by Paolasini et a/£3, is also consistent with the experi- 
mental findings. However we feel that their argument was 
rather scant and we try to present here a more in depth 
analysis of the implications of X-ray resonant scattering 
data, deferring a more complete treatment to a future 
work. We shall see that the space magnetic symmetry 
group of the AFI ground state plays an essential role in 
the discussion of the structure factor, since it can dis- 
criminate among different theoretical models. 

In the anomalous diffraction at the V K-edge, the 
structure factor F is given by 



^ = E/"^ 



.iQ-R,z 



(7.1) 



where the sum is performed over the eight V-atoms of 
the unit cell, /„ is the resonant atomic scattering factor 
(RASF), Q the momentum transfer at the chosen given 
reflection and Rn the atomic position inside the cell. 

We follow the numbering of V atoms as given in Fig. 
m, and divide the eight atoms in the monoclinic unit cell 
into two groups with opposite orientation of the mag- 
netic moment, 4 (u,v,w), 5 {—u,—v,—w), 2 (1/2-1- 
u, —V, w), 1 (1/2 — u, V, —w), and 8, 6, 3, 7, with coordi- 
nates obtained by adding the vector (1/2, 1/2, 1/2) to the 
flrst group. Here u = 0.3438, v = 0.0008, w ^ 0.2991 are 
the fractional coordinates of the atoms in unit of the mon- 
oclinic axig3 and the origin is taken between the atoms V4 
and V5 in the same basal plane. It is important to note 
for the following discussion that the be monoclinic unit 
cell (space group /2/a) is not primitive, so that, neglect- 
ing the magnetic moments, the two groups of atoms with 
their oxygen environment are translationally equivalent. 

At the (/i, k, I) monoclinic reflection, putting for 
brevity b = e^^(''+'=+'\ and_^ = 2Tr{hu ± hv + Iw), the 
structure factor F in Eq. (7.1) can be written as 



F 



(/8 + 6/4)6^* + (/6 + 6/5)6- 



-i(j>+ 



(7.2) 



H~iy 



{h + bh)e'*~ + (/i + 6/7)6- 
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The RASF /„'s in this equation-are tensorial quantities 
given by the following expressionEj: 



fnii^) 



^ nuj + Els -E^- iTn/2 



(7.3) 



The transition operator O is the electronic part of 
e-f(l+ik-f) with the usual notation for the photon energy 
u), the polarization e and the propagation vector fc, so 
that fni^) may bear two, three or four Cartesian indices 
according to dipole-dipole (DD), dipole-quadrupole (DQ) 
and quadrupole-quadrupole (QQ) transitions. Then the 
atomic scattering amplitude is obtained by saturating 
with the appropriate e'and k components. 

We can now use the operations of the magnetic sym- 
metry group of the crystal to transform the /„'s into one 
another. In doing so we remember that a symmetry op- 
eration changes the site index according to Table I of 
Section IV and simultaneously acts on the tensor com- 
ponents. This action reduces to a multiplicative factor 
±1 since P = C^ = a^ = f^ = E. Thus //„ = /„ for 
DD, QQ and //„ = — /„ for DQ transitions. Moreover, 
we find it expedient to split the RASF tensor into a sym- 
metric (charge) part and an antisymmetric (magnetic) 
parti: /„ = /+ + f- . In this way, f/± = ±/±. No- 
tice that in establishing the site correspondence and the 
transformation properties of the tensor components it is 
essential that the Wannier functions of adjacent hexago- 
nal planes be defined as in CNRla (Eqs. (2.4) and (2.7)), 
i.e., they should be related by the same symmetry oper- 
ation connecting the sites. For example, since the sites 2 
and 4 are related by the glide plane (Tf,, the definition of 
the Wannier function at site 2 should be obtained from 
that at site 4 by changing y into — y, if we take the ref- 
erence frame as depicted in Fig. yj 

Consider now the magnetic group C2h®T, which is ap- 
propriate to the ferm-orbital AFI phase of V2O3 as pro- 
posed by Mila et alH Defining A{(j)^) = (e*"** -h/e-*"^*) 
and using the symmetry operation / to connect sites 5 
and 6 to 4 and 8, and sites 1 and 7 to 2 and 3, and the 
operation Tai, to connect sites 2 and 3 to 4 and 8, we can 
express the structure factor as: 

Fhki = [A(0+) + A{r){-lffa,]{f+ + bf+) 

+ [A(0+) + A(r )(-l)''rab](/s- + bU) (7-4) 

Using now the peculiar fact that in V2O3 w w 0, defin- 
ing (p = 2Tr{hu+lw) and A{(j)^) = A{(p) and remembering 
that Tf^ = ifn, we finally obtain: 

Fhm = (e*^ + /e-^^)[l + {-lta,]if+ + bf+) 

+ (e^0 + /e-»0)[l _ {-l)>^a,]if, + bf^) (7.5) 

In this expression we assume that the operators have 
already acted on the site indices, so that they should act 
only on the tensor components. We notice that, consider- 
ing the /„ as scalars, we recover the crystallographic ex- 
tinction rules (the reflections with h + k + l odd (6 = —1) 



are lattice forbidden) and the selection rules for non res- 
onant magnetic scattering (for 6 = — 1, those with h even 
are magnetically allowed, those with h odd are forbid- 
den). The classification by Paolasini et alS3 of the lattice 
forbidden reflections into purely orbital and purely mag- 
netic was actually based on the non resonant magnetic 
scattering and might therefore be only approximate. 

Returning to the resonant regime at lattice forbidden 
reflections with h odd, one can easily see from Eq. ( |7.5| ) 
that there cannot be any charge (and therefore orbital) 
signal if sites 8 and 4 are in the same orbital state, since 
they are completely equivalent from tiie point of view of 
the lattice. Therefore what Mila et aliM are really saying 
is that the (1,1,1) reflection observed by Paolasini et al. 
can only be of purely magnetic origin, which however has 
nothing to do with orbital ordering. 

A similar argument is also valid for the C2®T symme- 
try of the RS(AO) phase. In this case, using T and C2 to 
relate the various sites, as shown in Table I, we obtain: 



Fhki = (e** 



icj) 



+ (e 



+ b{-lfC2e-'"^){b + 1)(/+ + b{-lff+) 
- b{~lfC2e-^'^){b - l)(/3- + b{-lff^) 

(7.6) 



and again it is impossible to observe a charge scattering 
when b = —1. 

We believe however that the claim made by Paolasini 
et al. for evidence of orbital ordering has some substance, 
despite their approximate classification. Indeed it is well 
known that the magnetic part of the tensor for DD tran- 
sition is only active in the cr— tt channel and its imaginary 
part can be shown to be proportional to the magnetic cir- 
cular dichroism in absorption. At the K-edge this latter 
is proportional to (L) , the averag&^f the orbital momen- 
tum operator in the final atate.E3'EZl Due to the fact that 
in V2O3 the moment liescl on the glide plane, a refiec- 
tion with respect to this plane change the sign of (i), 
so that af~{DD) = —f~{DD). This implies, from Eq. 
(7.5), that the DD tensor is magnetically active only at 



refiections with h even, but cannot contribute at Ji.pdd. 
The experimental evidence from Paolasini et aZ.E3'E3 is 
in keeping with this conclusion except for the observed 
intensity at 5464 eV, corresponding to Is — *■ 3d transi- 
tions, where at both types of reflections {h even and h 
odd) signals of comparable magnitude were detected. 

Notice that we expect a strong contribution from the 
DD transitions at this energy for h even, since the V- 
atoms are off-center with respect to the oxygen octahe- 
dra, which are furthermore distorted. At the same time, 
the ratio of the transition matrix elements between DD 
and DQ transitions in the case of V2O3 is estimated by 
Fabrizio et a/.O to be of the order of 7, in goodpSgreement 
with the experimental data by Goulon et alzA. There- 
fore the observed comparable intensities at 5464 eV for 
reflections with h even and odd cannot have the same 
magnetic origin. This argument also iavalidates the ex- 
planation proposed by Lovesey et als^ that the (1,1,1) 
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signal is due to the octupolar component of the QQ tran- 
sition, a tensor of rank 3, apart from problems related to 
its use of an approximate symmetry. At this reflection 
the signal should indeed be much smaller than at an al- 
lowed DD magnetic reflection [h even). Therefore much 
more work is needed to establish the exact implications of 
the resonant diffraction experiment on the nature of the 
AFI ground state of V2O3, since this complex question 
cannot be dismissed by too simple arguments. 

If we now consider the reduced group C2h{Cs), using 
TI and cTb to relate the various sites, again according to 
Table I, we find for the scattering factor : 

Fhki = (e'^ + bfie-'^)[l + b{-lfab]{h + bh) 
= (e^^ + 6/e-^*)[l + b{-lfa,]U^ + bf+) 
+ (e'0 _ bie-^^)[l + 6(-l)'^a,](/g- + bf^) (7.7) 

At reflections forbidden by the lattice (b — — 1) , the 



charge scattering behavior of expressions (7.5), (Z.E) and 
( [7.7I) is substantially different. As already emphasized in 
the case of the two groups C2/1 (81 T and C2 ® T there 
cannot be charge scattering, while for C2h{Cs), fs is dif- 
ferent from /4 , so that a charge signal might appear. 
In particular this is the case for the orbitally ordered 
state RS(ME), due to the different molecular states of 
the molecules straddling atoms 8-7 and 4-1. This fact 
would readily explain the comparable intensity observed 
at reflections with h even and odd, since the transition 
is always dipolar in nature, and stronger than that due 
to DQ or QQ processes. Again more work is needed to 
establish whether the predictions of the ground state pro- 
posed here for V2O3 are in keeping with the experimental 
flndings. 

Finally we want to point out that the RS(AO) con- 
flguration, due to the orbital ordering, can give rise to 
magnetic reflections of dipolar character, with 6 = — 1 
and h even or odd, as apparent from Eq. (7.6), since 



atoms 3 and 4 belong to vertical pairs in different molec- 
ular states. However this phase cannot be an eligible 
candidate for being the ground state of V2O3, due to the 
presence of T in the associated space magnetic group, 
that forbids magnetoelectricity. 



E. Conclusions. 

In this paper we have derived the strong coupling limit 
of an Hubbard Hamiltonian with three degenerate i2g 
states containing two electrons coupled to spin S = 1, 
and have re-examined the low-temperature ground-state 
properties of V2O3. Given the present experimental ev- 
idence in favor of a spin-1 state of the V'^^-ions, thejta 
is no doubt that the original suggestion byJ. Allen,EHl 
later taken up and developed by Mila et al&B, that the 
vertical molecules constitute the building blocks of this 
compound in all three phases of the phase diagram, is 
the key concept to understand the physics of the phase 



transitions observed in V2O3. In fact the molecular unit 
reconciles the existence of orbital degeneracy, required 
to explain experimental evidence ranging from inelastic 
neutron scattering data to the relaxation time of the nu- 
clear spins in NMR studies, with the spin-1 state of the 
V'^^-ions. Focusing on the AFI phase and analyzing the 
parameter space of the problem, we have indeed found 
that one has to consider two regimes, depending on the 
relative size of the intramolecular correlation energy and 
the in-plane exchange energy. These two regimes dic- 
tate the form of the variational wave function. In both 
cases we found minimizing solutions in a reduced region 
of the Hamiltonian parameters (hopping integrals versus 
-^) with the real spin conflguration. The molecular so- 
lutions are more suited to explain certain experimental 
facts, like, for example, the lack of detectable anisotropy 
in various spectroscopies in the paramagnetic phases of 
V2O3. However there remain two orders of problems 
which will require future investigation. 

The flrst one is the rather small stability region in 
the phase space of the parameters — versus jj- and the 
small stabilization energy of the RS phases compared 
to neighboring competing configurations (in the range 
1.5 -^ 2 meV, using the standard set of parameters). 
The AF(AO) phase, with its continuum degeneracy and 
the in-plane spin-orbital antiferromagnetic coupling with 
trigonal symmetry, is very reminiscent of the paramag- 
netic insulating (PI) phase at higher tempierature, setting 
aside the short range of the correlations^ Therefore one 
would expect an energy gap at least of the order of the 
transition temperature (15 meV), and we knowQ that the 
estimated Neel temperature of the AFI phase is about 
20 meV. The problem shows up again if one tries to cal- 
culate the exchange integrals for spin wave excitations. 
Even though, looking at Fig. |^(b), one can find a value 
of Hund parameter J such that the exchange integrals 
along 61 is ferromagnetic and about half the one antifer- 
romagnetic along 62 or S^^s found by inelastic neutron 
scattering measurements,LJ their actual values are off by 
a factor of five. Moreover it does not seem that longer 
range hopping parameters, neglected in our calculations, 

might cure this drawback. Since the energy scale is set 

2 
by JJ-, in order to get the energetics correct one should 

increase r by more than a factor of two. However this 
would be an easy way to escape a problem that has more 
profound roots. We tend to believe in fact that the diffi- 
culty steins from the mean fleld treatment of this highly 
correlated electron system on one side and on the other 
side from the neglect of the lattice degrees of freedom. 

The second order of problems regards the consistency 
of the RS phases with the experimental flndings of x-ray 
synchrotron radiation spectroscopies, namely the anoma- 
lous diffraction and non-reciprocal x-ray gyrotropy. We 
have argued that the comparable intensity of the ob- 
served signal at the Vanadium K-edge in resonant diffrac- 
tion at both reflections with h even and odd cannot be 
interpreted as being purely magnetic in origin, but finds 
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a natural explanation if one assume that the ground 
state magnetic symmetry group is such as to allow charge 
scattering, i.e., if there is the appropriate orbital order- 
ing. This group is the same as the one indicated by the 
non-reciprocal spectroscopy. However none of the sta- 
ble RS phases found admit this group, so that they only 
allow resonant magnetic scattering. This is the reason 
that prompted us to indicate the excited configurations 
RS(ME) as possible candidates for the correct ground 
state of the AFI phase, invoking an optimal coupling 
with the lattice degrees of freedom. This assignment 
makes also a definite prediction for a transverse magneto- 



electric effect in this phase which could be subjected to 
experimental test. 
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APPENDIX A: THE EFFECTIVE HAMILTONIAN IN TERMS OF FERMION OPERATORS. 



Here we report the results of the calculations performed to derive the effective Hamiltonian in terms of the fermion 
operators. 



The explicit form of the projection operators X^^ , {i = 1, 2, 3) defined in Eqs. (2_^), (2^) is given by: 
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For simplicity we omitted the site index as all operators act on site j. 

Dividing iJoff into three parts, according to the X^^' definitions, we obtain the following expressions: 
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3{U2 + 2 J) ^ , 

^ ' ijcT nn' 



/ J / J ,/_v V ^mni'Jl,-'- ^?nm"Jl,-l- ^m'va") 



27 ^72 \^infj'--in' a "'jm' (7 ""jm" a ^^ c\ ina^iT^ '^ J^^ <^ 3''^ '^ ) 



.nm.m'n' f + „. , „+ „. „. „ _ i p+ r- / - r+ r- T? ■ // + r+ r / r+ r- -r? ■ // "I 
y :/« y-ina'^^ri' iy<^jjn' a^jma ii-jm" a ^ '-ina'^in' a<^jmt g<^jm(j ii'jm" a T t-^^gL-m' a<^jm' a^jma "-jm" a ) 



"^(TT -^9 7"'! / ' / ' / y V "nim'jl,-'- "mni"j(,-'- "m'm"; 

ijcr nji' mm'm" 



.nm.rnn' f+ + + 



'^^in<j'^in''y''^jm'aCj^,,gCjm"a + C.ins'^in' aCj^^i „Cjm' a^jr 






+ ^--ina^iri' aCj^i gCjmanjrn" a + C;^,^^Cini aC^^-^i ^CjmaCjm" <jCjm" a + C;^,^^Cini aC^^-^i ^CjmaCjm" aCjm" 



' '^ina^i'n'(jCjjj^,^Cjm(7Cjm"(7Cjm"a + '^ina'^in' (jCj^i ^Cjmanjm" a ) 



(A5) 



H-(3) _ 



^EEEd-^"-) 



2(t/2 + 4J) 



2J(T nn' mm' 



j.nm.mn'f+ f,. ,„._„., i f,+ ^. ,_„ . „. 
y J« \'^in(y^^ri'aii']ma"']m'a ^ '^ing'-in'ai'-jmaii'jr 



'^inrr^in.' aCjj^gCjmanjm' a C^j^^Cin' aCjj^^Cjmanjm' a 



9(71 I 1 j\ / , / , / y V-'- 0mm' j(, J- dmm"j(,-L dm'm"J 
ij(T nn' mm'm" 



.nm.m n { + 



fnm.mn/ + r- , c'^ C ,-n ,i - r+ C / - r+ C , V n 



^ina'^in' trC-jma'^jm' crl^jm" a + ^^inS^in' o^jmcr^jm' a'^jm" cr) 



(A6) 
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APPENDIX B: COMMUTATION RELATIONS. 

In the following wc evaluate the commutators used to derive the effective Hamiltonian in section 11. 
As usual, given two generic operators A and B, we define [A, B] = AB — BA and {A, B} = AB -\- BA 



'"ma ""m a ""ma i 



^ra'" a"'\ "^ ^mm'" ^aa'" ^ma^"^^' '^''^'m" a" 
y^m'm'" ^a' a"''^ma^jYi> (j''^m" a" "^~ ^m" m'" ^a" a"''^ma'^7n' a' ^^ 



+ 

'm" a" 



(Bl) 



\Cm"' 

-\-Sm' 



' ; {""ma ""tn' a' ""m" a i 



^m'" a"'\j — ^mm"" ^aa"" \Pmm"' ^aa'"'^' 



m'a' ""m," a" 



ma ""m,' a' 



' ^a' a'" C-maC-^^/TljYi" (jii -|- Ofn" m'" ^a" a'" ^\ 

iin I nil \ n ii inn n iiiT) C i if —I— n i nin i iiiTI T) n n 

"U(j'(j"" yJiri"'m"' ^a" a'" "^ma^'m' a' '-^^^^'f (j" \^ ^m'ln'" ^a' a'" "ina""m" a" 
0(j(j"' C^j^^CjYi'fjiTlm" a" ) ~r ^'m"ra""^a" a"" \Pm" m'" ^a" a"''^ma'^m' a' 
'^a'a'"^maC / /C'(ji"a" ^mm'" ^aa'" ^rn 



(B2) 



jYia 'f^ (^ ^Tn" a" 



[nm" a" C^^Cjn' aC^^i ^Cjna ■, C^n>^ni\ — 0^" m'" Oa" a'" C^n ^ffC^^jCjn' aC^-^, ^Cjna 
'i^TijYi" (j" C^^^^C,^^! ^lyOj-fii jjiiii 0(j(j"' CjYi^^ " 0mm'" ^aa'" (^m' a ) 



(B3) 



\Cra"" a"" -J \P"m" a" ^rna^ma^fn'<7^fn' (7 t ^rn'" a"'\i — ^m" m"" ^a" a"" \Om"m"' ^a" a'" (^rria^ma^'' 
-rOfn' m'" C-m" a" ^ma^ma \Pa"' aC-m'a 



m' a^m' a 



fCm'a)) 1 0mm"" ^aa"" \0' 
' Oa"a"' ^m,"a" ^ma^m'aC' 
-Om'm"'^m"a" ^max o'a'" ^m'a ^t 






a'^-mi < 
m"m" 
m'm" 



I'm" a" '-ma 



cr\Oaa'' 



m mi 

m'a) ~r Oinmi"" Oaa"" \0', 
aa'" C^m'a ) ) 



s. 



aa ^m a ) 



m" m"' 0(j" a"' ^^mi" a"^ma^'rfi' a^mi' a 



(B4) 



i^m" a" ^ma^mia^Tn' a^^m' a ^ ^>iji"' (j"'\ ~^ '^rn"m'"^(T"cj 
'v'^m" a" (^rria^mia^Tn'mi"' \0(j(j"' Cjjii ^ Oaa"' ^^m' a) 



7" ^mia^ma^fn' a^m' a 



(B5) 



\Cm""a"" J \P"m" a" ^m(j^m' a^^ ^^ma-, C^tt ^ni\ j — O^a" mi"" Oa" a"" \Pm"m"' Oa" a"' ^mia^"^^' <^^mi' a^"^ 
~rCm" a" ^rna^mi' a\ '^^''^^" '^'^"' ^'^^^ Omim"' Oa"' a^^m' a ) ) 1 Ojnm"" 0(ja"" y^mi" a" C-^ ^ 

X yOj-fi'jYi'" 0(j(y" 



-'mm '-^aa ^-711 a ) 



Om'm"" 0(7(7"" \0 



"'mi mi ^'a a rri" a" rn <j'-jyif Q-'-ma ) 



1/" 0(7" (7'" 



m" mi'" ^a" a'" L-^^//^//(--^^(j-(--rn'cr*--rn,(T 



nC' 



1 '^mi" a" ^ 'Y-[i(7 \0 jYi' j-fiiii 0(7 (7 III CjYi(7 0mm"' Oaa"' ^m' a) j 



(B6) 



APPENDIX C: THE EFFECTIVE SPIN-ORBITAL HAMILTONIAN. 

As anticipated in section III, the spin orbital representation of the iiTeff is quite cumbersome, so it was not reported 
in the main text. Nonetheless we think it could be useful to write down its complete form because, despite its apparent 
complications, it becomes very simple to handle for real materials, where, due to symmetry considerations, many of 
the terms could be zero. 

Starting from the expression of Appendix A, we introduce the spin-pseudospin representation defined in section III 
and write down H^s in the following form: 



H. 



(1) _ 1 1 V^ V^ + 

ij nn' 






\/2 -^ ^■'^ 2 ^^ ^* ^ -^ ^/2 ^^ ^^^ /2 -^ -^ /2 -^ -^ 



[S^ ■ Sj + 2] (CI) 



H^'J = - 



1 1 
UU2 + 2J 



ij nn' 



f^t^'(r!z + T,. + 6) + Itftfii. - r,, + 6) + ilfi^'ii - 4) 
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^^ (-^^i^^^ ^^^^ 






(C2) 



,«2,3n' jLr+T- - x/2r- T- 



^^ V \/2 -^ -^ ■^ 



•'ij "ji 



t"ft?r' 4r,.rr-V2- 



V2 



J2'j 



'j 'jz 



[l-S,- s,] 



cn 



1 1 



4 C/2 + 4 J 



ij nn' 



^in^i'' 



lt^lt]f {I ~ T,.){2 + r,.) + It^tf {2- r,.){l + T,.) + t^tf'Tf^ 



}_fnl,2n' - - L/"l/3n'_"-7-- + i/"2.1n' + + L+"2,3n' + 



(C3) 



j.n3-/-ln' ^+ 



j-n3j-2n' 






[1 - 5. • 5,] 
The summation on n and n' must be taken according to the rules of the fohowing table: 



n-1, 




= 1 


n-1, 




= 2 


n-1, 




= 3 


n = 2, 




= 1 


n — 2, 




^2 


n = 2, 




= 3 


n = 3, 




= 1 


n = 3, 




= 2 


n = 3, 




= 3 



CitiC,„/ = i(l-r,J(2 + riJ 



'"'in^in' — 2'^i '''i 



V2^ ^^^ 



C^'-n Q^ 



CitiC„i' = 5(2-Ti^)(l + ri^) 



^in'^i'' 



'^in'^i' 



^2'^^ '^'''■ 



^T- T+ 



"^' 



-'in'^in' — ' iz 



If we perform the summation over n and n' explicitly, we get: 



r(i) 



1 1 



<=«■ " 3 t/o - J 



— 5][S^..5,+2] 



(41)2(2 -T.,-ri^)r,.(r,. + l) 



{i^f{2 + r.z-rf^r,z{r,z-V) + ^(1^frt^(l-r^^ 



1.11 



4/^^l(rr^ + ^^)(i 



r,, -t2^) + — tlltlf(T-T,,+T,,T+)(l-X- ~2 



^/2 



jz '^j^y 



o*'/*??^ ^ 44 + ^*'^?' (^ ^ ^-^"' + ^' ^' ^^ ^- 



2 y J* * « J J /o 



'ij '•ji V'i 'i 'iz'j ^ 'i 'i 'j ijz) 



aU33 t - H 

""ij ^]i \'i '*z 'jz 'j 



') + \{t\^f[^rl+r. 



22 /-----I + + + + 

^zT-jz - T^zTjz - 2 V» ^» '^j ^J + '7 ^J V V 
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+ :j^^ irir-rjr,. + ^^r^.r^ + (^.^ + ^^.)(1 - r,. - r^)) 



1 

V2^ 



+ —fK^f3^f^ k^i ^^ ^i^'^J- + '^t'^t'^t'^-J'- + (^»2^i^ + ^» ^*2)(1 + ^J^ ^ '^1^^) 



ij ji Vi '''i ^" ''"i ''"i j V- ^ '^3^- '^jz) ^ij'^ji Vi TizTjzTj + TizT^ Tj Tjz) 



+ 7;it\f? ((2 - 34 - nz + 2nzTl + 2^4) - ij-nzT-T,, + T,,r+T,,T+)) 



ilf4' ( r^nzT+Tjz + nzT+T,zTT- + i(r+r+ + r, r, )(1 - 2t^ 



•]z<j ^2* 



jz- 



+ 7f *i^j'^' (^^ ""* ^/^^^ + ^^^t^l^-^3 



-^Tf^'^f (^-^^+^rn.)(l-2r_ 



^i?^,'f(r..rr+r+r,,)(l 



^j-z - r%) - -^^ {nzT- + T+nz) (1 - 2- 



jz^ 



tjpf, {nzT^ T+T,z) + -{tf^f ((2 - ZtI + T,, - 2nzT% + 2tIt%) - {t+T,zT+Tjz + r,,T^. T,,T, )) 



H. 



1 1 



(2) 

""^ 12 C/2 + 2 J 



E[i-^'-^.] 



(i-/)"^ ( 3 - Ti^ - rf^ - -TiyTjz{l + 2Tj2 + TizTjz) 



+ mY[^ + nz-rt- 



\t,zTjz{\ - 2nz + nzTjM + (t33)2 (4^2^ „ rljl) 



\t\}t]l {r-r- + T+T+) (2 + T,z + T%) - ^t"t]? {{T-r,z + t,,t+) (2 + r,, + r/j) 



it"t22 (' 



^^^0^3^ r. r^T^ T^ + T^ T^ T, T, + 2(t, T, - T, T, ) (t^- T^- - T^- T^T 



-)) 



^q {rrr,z{2r+r+ ?,rjr-) + T,.T^(2r-r- - 3r+r+)) 



+ :^i!^?' (TrT,,(2T,,r+ - irrr,,) + t,,t+(2 



Tj -Tjz - 3rj^T+)) 



ol^'i")' [(^^^^(^/^/ - 3^7 rr) + rrrrr-rr) + (12 - ^rf^ - 4^ + n^r^z)] 



^4Mf (t-T,,(2t-T- - 3T+T+) + T,,T+(2t+T+ - 3T-T-) + {nzT- + T+T,,)(2 + T,, + T^j) 
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+ 2^1!^^^ ^^^^^ + ^^^^)(2 - r,, + rjj) 



- ^ii^'f {n,Tr{2T-T- - 3T+T+) + t+t,42t+t+ - 3Trrr) + {n,r+ + T^n,){2 - r,, + r^)) 



+ t\lAl {nzrr{2TrT,, - 3t,,t+) + r+T,,(2T,,T+ - 3rrr,,)) 



JZ'j ■J'j <]Z 



(tljf ((4 - 2r,^ + Tj^rj*^ + r^^T^J + r^ ^^(t^. t^-^ - 3tj^t+) + Ti,T+Tj,T^) 



V2 



i^fi?? (n.rr(2r+r+ - 3x7^-) + T+T,^(2r-r- 3r+T+)) 



+ tlfl^f (r,,Tr(2r,,T+ - 3t-t,,) + r+r,,(2r-r,, - 3t,,t+) + {t+t+ + TrTr){2 - r^J) 
- ^tlft'f ((^-^^ + ^r^-)(2 - r^J) - ^i'M ((^..rr + r+T,,)(2 - r,, + r^J) 
^4'*f (r,.Tr(2r+r,, - 3t,,t-) + r+r,42T,-,Tr - 3t+t,,)) - ^2^2/^33 ((r.^rr + t+t,,)(2 - r^J) ; 



+ (4')' ((4 + 2T^. - n.rf, + tIt%) + T+r,,(T+Tj, - Sr^.T^- ) + r,,r, t,,t^ ) 



H. 



(3) 



1 1 



off 



4 t/2 + 4 J ^^ 



[1 - 5. • S',] 



\^'^^ 



{t)\f(2-r,^-Tl)(2-T,^-rl) 



^f{2 + n.-rl){2 + r,.-rl) + {i^f {rgrf^ 



+ 7^t\}t]l{riTr +T+T+){2-T, 



rl) 4t"t!?(7 



]z 'jz) /K'^ii'^ji^'i ''iz^'izli )y^ '> jz 'jz) 



'n^ + n^T+)(2 - Tj^ -tI 



+ 4%*J« Vi ^i '''j '^J +"^5 'J 



+ + + +^ _ J_+ii+23 (' - - - , + + +^ 



4'if {r, n,T^ Tjz + Tj,T+n,T+) + ^{tlff ((2 - n, - tI){2 + Tj, - T%) + T, r, r; 



n ij ji Vi "^^ ' j '3^ ~^ ' 3^ ' j '*^'i / ^ o V-Jj 



V2 



*J,^*i^ ('^i ^^ ^J^'^7^ + ^i^^i^^7 ^^'^ + (2 - Tiz - Tfj(T+Tj2 + Tj^T. )) 



-ti2t22(^-^- + r+r+)(2 + T,, - r^,) 
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- ^^^f^^i^ (^« ^' ^^^^- + ^^^^^1-^j + (2 + rjz - t^J{t^ Tj, + Tj,T+)) 



■ ttft^f (r, TizTTTi, + Ti^T^Ti^T^ ) + {t\^f ((2 - T,^ - T^ jrj^^ + T,- t,:^t,-^t; 



il2.33 



<i 'tzijz'j 






T- T- TiVT- + T- r- r- Ti 



)+43if(T^ r,zr,,r, + t,,t+t+t,, + tI{t^ r^ +r+T+)) 



^ y J» *> i 'iz^'iz'i )'jz /nVJi^i 'iz^'tzli )(^^ljz Ijz) 



V2t^^tf,{T+nz + T,zT-)t, 



1 



O ij 3^ Vi ''''^^■'''j ''jz + 'T^zTj- TizT^ 



For the sake of brevity, we have not symmetrized the 
above expressions over i and j, so that we wrote, for 
example, the expression ti^Tj*" + T2zT^ + t^ 



i-lz 



as 



E 

U = l,2 



nzTj +r,z 



instead of 



ij=l,2 



Note that even if in principle weglapuld have 81 terms 
(3''), we have assumed, as in CNrJ3'EI that V^^' = i^'™ 
so that some regrouping of the terms has been possible 
(each of the three if^g contains only 21 terms). Unfor- 
tunately, as anticipated in the text, the fact that there is 
no conservation law for the pseudospin quantum number 
Tzi during each hopping process, prevents from having 
more symmetrical expressions for the orbital part. 



APPENDIX IV: ACTION OF Heff ON 
MOLECULAR AND ATOMIC STATES. 



In this appendix we show how to derive Eqs. ( p.5D , 
( |6.7[ ). For a given horizontal bond (for example, 8\ in 
Fig. l2|) nine two-site orbital states are involved in the 
product state (3.4). They are hsted below: 



-rl)rl^ 


Ttr^zT^zT^ ) 


II l>ac = 


-l)a|-l>c 


II 2)ac = 


-l)a|0)e 


II 3)ac = 


-l)a|l)c 


II 4)ac = 


0)a|-l)c 


II 5)ac = 


0)a|0). 


II 6)ac = 


0)a|l)c 


II 7)ac - 


l>a| - l)c 


II 8)ac = 


l)a|0)c 


II 9)ac = 


l>a|l)c 



(4.1) 



If the spin configuration of the bond is ferromagnetic, 
then after some simple algebra we find: 



K 


l>ac 


K 


\2)ac 


K 


\i)ac 


K 


4)ac 


K 


5)ac 


K 


6)ac 


K 


7)ac 


K 


8)„e 


K 


|9)ae 



2t^+2x^ 
U2-J 



l)ac 



/3 +g +x +e 

U2-J 
U2-J 

U2-J 



^ U2- 

^/ ac 
4)ac 



,7 II ^)ac 
U2-J I 



7). 



2t''+26I^ I 
U2-J I 



'^1 ac 

±xL\ 



(4.2) 



(72-./ 
U2-J 



7)ac 



2a/3 
U2-J 



II 3). 



ct+a^+T^ 



U2-J I 

.2x!±2»! 1] q\ 



/ ac 
2X^ 

U2-J 



1). 



whereas for the antiferromagnetic configuration we ob- 
tain: 
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Hf, II l)ac 



-.2 I ^2 I n/]2 



'5i 



^.1 



^.1 



>> (C/2+4.7) 



3a^+3cr^+4r^+4x^+66>^ ^ ii -, v 

3(L'2+2J) / II -^/"c 



+ ( 



a/3 



a/3 



Hf, II 2),, = -( 



(72+4,7 C/2+2J 



2 I _2 I ,,2 I n2 



) II 9)ac 



_ _^af+Tf+X_+£ 
((72+47) 



3((72+27) 



^.1 II 3).c = -( 



2 I —.2 I 2 I /i2 



_ ^ g-^+r^+X^+e 

((72+47) 



2 I 0/^2 I Q_2 I f-_2 I Q,,2 I c/32 , 



2a"'+2/3''+3CT''+5T^+3x +56 
3((72+2J) 

»/3 II 7)a, 



2). 



3>c 



3((72+27) 



^.1 II 4)ac 



^ a^+T^+X^+e^ 
^ (C/2+4J) 

+2ct' 

3((72+2J) 



4). 



rr^ II c;\ - _/a^+^f+2x^ 

-H^. II O^ac - (, ((72+47) 



3a^+3j3^+4r^+6x +46>^ n, n r\ 

3((72+27) / II "^/ac 



6),, = _(^!±i!±x!± 



rT-A II y\ _ / g- +-^ +X ^ 

-"5i II 'Mc — I ((72+47) 



(C/2+4J) 

2a''+3/3^+2g^+5T^+5x^+3e'' 
3((72+2J) 

2 I -.2 I , 2 I /i2 



_ /£f+r++xl+e 



2a^+2j3"'+3CT"'+5r^+3x +56"' 
3((72+2J) 

2a/3 



6). 



7). 



+ 3((72+27) I! 3)ac 



^2 I _2 I „2 I /]2 



_ _//3f+rj+2(f+ej 
~ V ((72+4J) 



2a^+3/3^+2g^+5r^+5x +36 
3((72+2J) 



^.1 II 9).c = -( 



g^+g-^+2T'' 
((72+47) 



3/3^+3(T^+6T^+4x^+4g^ N, II Qx 

3((72+2J) / II -^/ac 



a/3 



af3 



(4.3) 



(4.4) 



^(72+47 (72+2J-' II -'-/ac 

As in the text, we have retained the hopping integrals 
6 = t^^ and x = ^^^ which are zero in the corundum 
phase and can be different from zero in the monochnic 
one. 

Given the form of the molecular wave function in Eq. 



(5.5), these are the only terms which survive when taking 
an average of the kind shown in Eq. (6.4). 
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